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ABSTRACT 

Saltwater  marshes  and  wetlands  are  important  buffers  at  the  land-sea  interface.  Among 
the  most  biologically  active  ecosystems  on  Earth,  natural  and  man-made  wetlands  are 
important  interceptors  of  pollutants  and  nutrients  bound  for  the  coastal  ocean.  The 
transport,  dilution,  and  deposition  processes  occurring  within  the  marsh  are  key  factors  in 
determining  this  interception,  and  these  are  in  turn  determined  largely  by  tidally  driven 
flows  as  influenced  by  marsh  vegetation  and  other  physical  characteristics.  Vegetation 
type  and  density  are  of  primary  importance  in  these  processes,  both  for  pollutant/nutrient 
uptake  concerns  and  in  determining  hydrodynamic  characteristics  of  the  marsh. 

This  study  examines  the  effect  of  vegetation  density  and  ambient  flow  on  diffusivity 
within  a  tidal  marsh  canopy,  specifically  Spartina  alterniflora.  Vegetation  densities  from 
0-1.4%  stem  coverage  and  flows  from  2-12  cm/sec  were  investigated  using  Rhodamine 
WT  tracer,  with  resultant  measured  dififusivities  ranging  from  approximately  0.5-3. 0 
cm^/sec.  Diffusivity  was  foimd  to  be  a  strong  function  of  ambient  current,  but  a  much 
weaker  function  of  vegetation  density.  Presence  of  vegetation  caused  transverse  and 
vertical  difflisivities  to  be  essentially  isotropic  over  all  vegetation  densities,  as  opposed  to 
the  order  of  magnitude  difference  found  in  earlier  non-vegetated  studies.  Only  slight 
vegetation  coverage  was  foxmd  to  be  necessary  to  produce  this  isotropy,  with  little 
additional  change  as  stem  density  increased. 
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1.  Introduction  and  Motivation 


Pollution  of  the  coastal  ocean  is  a  problem  attracting  increasing  attention  as  coastal  zones 
continue  their  rapid  development.  This  development  directly  impacts  coastal  water 
quality,  and  hence  coastal  ecosystems.  Few  would  disagree  that  maintaining  coastal 
water  quality  is  critical  to  wise  management  of  any  coastal  region,  although  there  is 
disagreement  on  the  magnitude  of  the  threat  and  the  degree  of  action  required  to  protect 
coastal  ecosystems  from  pollutant  damage.  Reducing  pollutant  inputs  to  coastal  waters  is 
one  way  to  minimize  this  damage,  and  both  natural  and  constructed  wetlands  have  been 
shown  to  be  very  effective  in  making  this  reduction.  This  chapter  provides  a  brief 
synopsis  of  the  pollution  threats  to  coastal  waters  and  the  role  that  fringing  wetlands  may 
play  in  mitigating  the  effects  of  these  pollutants  on  coastal  ecosystems.  A  clearer 
imderstanding  of  this  role,  along  with  the  potential  improvement  in  wetland  functioning 
that  may  come  with  this  improved  imderstanding,  is  the  ultimate  goal  of  this  research. 

1.1  Threats  to  the  Coastal  Ocean 

The  overconcentration  of  population  and  commerce  in  the  relatively  small  land  area 
bordering  the  ocean  presents  the  largest  threat  to  coastal  waters.  The  specific  components 
of  this  threat  may  vary  widely  from  locale  to  locale,  although  deteriorating  coastal  water 
quality  is  a  concern  common  to  all  areas.  This  deterioration  may  have  severe  impact, 
with  effects  ranging  from  closure  of  shellfish  beds  for  harvesting  to  closure  of  beaches 
due  to  medical  waste  disposal  or  toxic  algal  blooms.  Loss  of  habitat  in  ecologically 
productive  areas  is  also  a  concern,  as  is  increased  erosion  and  escalating  storm  damage 
costs  that  this  increased  coastal  development  brings.  The  Committee  on  Environment 
and  Natural  Resources  (CENR),  one  of  nine  committees  of  the  National  Science  and 
Technology  Coimcil  (NTSC),  summarizes  all  these  effects  in  their  list  of  ‘major  coastal 
environmental  issues’  (National  Research  Council  (NRC),  1994): 
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•  Eutrophication 

•  Habitat  modification 

•  Hydrologic  and  hydrodynamic  disruption 

•  Exploitation  of  resources 

•  Toxic  effects 

•  Introduction  of  nonindigenous  species 

•  Global  climate  change  and  variability 

•  Shoreline  erosion  and  hazardous  storms 

•  Pathogens  and  toxins  affecting  human  health 

Although  the  NTSC-CENR  did  not  prioritize  among  the  major  coastal  environmental 
issues  it  identified,  other  groups  have  attempted  to  rank  these  issues  in  terms  of  the 
threats  they  pose.  For  example,  the  Committee  on  Wastewater  Management  for  Coastal 
Urban  Areas,  a  committee  of  the  Water  Science  and  Technology  Board  (WSTB)  of  the 
NRC,  identified  nutrients,  toxic  organics,  and  pathogens  as  “high  priority”  on  a  national 
level  for  managing  coastal  water  quality  (NRC,  1993).  While  they  go  on  to  say  that 
different  constituents  may  pose  more  severe  threats  on  a  local  level,  these  three  are  most 
important  on  the  larger  scale. 

Although  toxic  organics  and  pathogens  may  be  significant  hriman  health  concerns, 
excessive  nutrient  addition  arguably  poses  the  greatest  threat  to  coastal  oceans. 
Eutrophication  is  one  of  the  primary  dangers  facing  estuarine  ecosystems  (Land  Margins 
Ecosystems  Research  (LMER)  Coordinating  Committee,  1992,  NRC,  1993,  NRC,  1994). 
Although  coastal  and  estuarine  waters  are  aheady  among  the  most  heavily  fertilized 
environments  in  the  world  (LMER  Coordinating  Committee,  1992),  primary  production 
in  these  areas  is  still  limited  by  available  nitrogen  supply  in  many  cases.  Primary 
production  is  closely  tied  to  all  other  ecological  processes  at  work  in  the  estuaries,  so 
maintflining  it  within  natural  bounds  is  key  to  preserving  the  ecosystem  as  a  whole.  A 
large  body  of  scientific  evidence  demonstrates  the  direct  relationship  between  nutrient 


inputs  and  primaiy  productivity,  with  nitrogen  being  the  nutrient  of  most  concern 
(Valiela,  et  al.,  1990). 

1.2  Role  of  Wetlands 

Although  wetlands  may  play  a  significant  role  in  reduction  of  all  pollutant  fluxes,  most 
data  available  concerns  their  nutrient  removal  capabilities.  For  simplicity,  this  topic  will 
be  emphasized  in  what  follows,  although  many  of  the  conclusions  reached  may  be  valid 
for  other  types  of  pollutant  inputs  as  well. 

Fresh  and  salt  water  wetlands  present  in  many  estuarine  regions  can  have  a  significant 
impact  on  the  flux  of  nutrients  entering  coastal  waters,  and  hence  on  the  estuary’s 
susceptibility  to  eutrophication.  Constructed  and  managed  wetlands  have  been  used 
extensively  for  natural  treatment  of  sewage  and  other  effluents,  and  the  natural  wetlands 
and  marshes  surrounding  an  estuary  serve  much  the  same  function.  Salt  marshes  are 
highly  complex  and  diverse  ecosystems,  with  biological  production  rates  hi^er  than 
those  found  even  in  intensively  used  agricultural  lands  (Teal  and  Teal,  1969). 
Additionally,  hinging  tidal  marshes  have  been  shown  to  more  than  double  soil  contact 
times  of  discharged  groundwater  relative  to  unvegetated  subtidal  sediments.  This 
increased  contact  time,  coupled  with  die  hi^  organic  carbon  content  of  wetland  soils 
(critical  for  denitrification),  translates  into  greater  opportunity  for  immobilization  of 
nutrients  prior  to  export  to  coastal  waters  (Harvey  and  Odum,  1990).  The  high 
productivity  and  soil  contact  times  of  these  marshes  make  them  very  effective 
interceptors  of  nutrients  bound  for  the  coastal  oceans. 

This  interception  occurs  through  a  number  of  mechanisms  which  vary  in  importance  from 
marsh  to  marsh.  These  processes  generally  impact  only  the  dissolved  inorganic  nitrogen 
(DIN),  the  dissolved  organic  nitrogen  (DON)  passes  throu^  the  marsh  essentially 
unaffected  (Valiela  and  Teal,  1979).  The  DON  is  less  energetically  useful  than  the  DIN 


to  organisms  within  both  the  marsh  and  the  coastal  ocean,  however,  so  effects  from  its 
input  are  less  severe.  Nitrogen  uptake  by  wetland  vegetation  is  generally  an  important 
sink  for  nitrogen,  and  nitrogen  fixation  by  bacteria  and  algae  also  make  significant 
contributions.  Data  from  constructed  wetlands  indicate  a  range  of  1 1-94%  nitrogen 
removal  capability  depending  on  the  degree  and  type  of  vegetation  present  (Reed,  et  al., 
1988).  Natural  marshes  exhibit  up  to  64%  nitrogen  removal  capability  (Valiela  and  Teal, 
1979),  indicating  that  some  enhancement  of  their  performance  may  be  possible. 


Much  of  the  nitrogen  not  intercepted  and  retained  within  the  marsh  ecosystem  is 
converted  in  form  before  export  to  the  coastal  ocean,  often  to  forms  which  reduce 
eutrophication  potential.  80-95%  of  the  DIN  entering  the  marsh  is  in  the  hi^y 
energetically  useful  nitrate  form,  but  the  majority  of  this  is  converted  to  particulate 
nitrogen  before  export  from  the  marsh  (V aliela  and  Teal,  1979).  This  particulate  nitrogen 
is  much  less  biologically  available  for  stimulating  primary  production  than  is  the 
dissolved  form  (Valiela  and  Teal,  1979),  so  eutrophic  effects  are  further  minimized. 

1.3  Description  of  Research 

The  preceding  discussion  delineates  the  critical  role  that  wetlands  may  play  in  reducing 
nutrient  (and  potentially  other  pollutant)  fluxes  to  the  coastal  oceans,  but  does  little  to 
quantify  specific  effects  of  various  wetland  characteristics.  This  is  very  difficult  to 
quantify  given  the  incredible  complexity  of  salt  marsh  ecosystems,  and  may  well  depend 
on  nutrient/pollutant  type  and  other  marsh-specific  attributes,  hr  general,  this 
exchange/conversion  process  is  governed  by  diffusion  from  marsh  sediments,  uptake  by 
plants,  tidal  flushing  characteristics,  and  other  factors  (Valiela,  et  al,  1990).  Vegetation 
has  a  major  impact  on  these  processes,  since  uptake  and  tidal  flushing  characteristics  are 
greatly  influenced  by  vegetation  type  and  density.  This  research  focuses  specifically  on 
the  effects  of  vegetation  density  on  diffusivity  within  the  marsh.  While  this  is  just  one 


parameter  in  a  combination  of  vegetation-related  parameters  influencing  the  marsh's 
nutrient  removal  effectiveness,  improved  understanding  of  this  effect  may  lead  to  better 
understanding  of  other  processes  occurring  within  the  marsh,  and  may  ultimately  lead  to 
ways  to  improve  wetland  performance. 


Equation  1-1  gives  a  simplified  eiqjression  for  the  transport  and  conversion  processes 
occurring  for  a  given  substance  within  the  salt  marsh.  All  uptake  and  conversion 
processes  have  been  combined  into  a  single  rate  constant  (.K)^  and  otherwise  the  substance 
is  assumed  to  behave  conservatively.  Advection  is  governed  by  the  3-D  currents 
throughout  the  marsh  (u,  v,  and  w),  while  diffusion  is  controlled  by  the  3-D  diffusivities 
(D»  Dy,  and  A). 


dC  dC  dC 
-f-  u  — r— +  V  -1— +  w 


dt 


dx 


dy 


dz 


d(  dC\ 

kc+—\d  — 
dxK  dcJ 
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dC^ 


ec'' 


Equation  1-1 


The  relative  importance  of  advective  versus  diffusive  transport  depends  on  the  scale  one 
is  interested  in.  At  larger  scales  (such  as  the  marsh  as  a  whole)  advection  dominates 
diffusion,  and  the  exchange  process  could  reasonably  be  thought  of  as  wholly  a  tidal 
advection  problem.  For  example,  to  determine  the  amoxmt  of  nitrogen  entering  Buzzards 
Bay  from  Great  Sippewissett  Marsh,  one  need  only  monitor  tidal  exchange  volumes  and 
nitrogen  content  of  exchanged  water  at  the  mouth  of  Sippewissett  Creek.  (This  marsh  is 
the  study  site  for  this  research,  a  more  detailed  description  is  found  in  Chapter  2.)  As  the 
scale  of  interest  is  reduced  diffusion  plays  a  larger  and  larger  role,  with  the  nitrogen 
balance  for  Great  Sippewissett  Marsh  as  a  whole  providing  an  excellent  example. 
Groundwater  flow  provides  roughly  90%  of  the  total  nitrogen  supply  to  the  marsh,  with 
most  of  this  flow  diffusing  up  through  the  bed  in  sandy  creek  channels  or  from  a  number 


of  springs  located  around  the  periphery  of  the  marsh  (Valiela,  et  al,  1990).  Groundwater 
is  unlikely  to  enter  the  marsh  in  significant  quantities  in  any  highly  vegetated  areas,  since 
these  areas  generally  have  a  much  lower  hydraulic  conductivity  (due  to  the  higher  organic 
content)  than  other  areas  of  the  marsh  (Valiela,  et  al,  1990).  Vegetated  areas  are  critical 
for  nitrogen  removal,  however,  with  over  95%  of  bacterial  nitrogen  fixation  and  roughly 
50%  of  the  denitrification  for  the  marsh  as  a  whole  occurring  there  (V  aliela,  et  al,  1990). 
These  processes  together  account  for  over  85%  of  the  marsh's  nitrogen  removal  capability 
(Valiela  and  Teal,  1979).  It  is  therefore  necessary  to  accoimt  for  the  transport  of  nitrogen 
firom  the  non-vegetated  regions  where  it  is  introduced  to  the  vegetated  regions  where  it  is 
utilized.  This  transport  is  both  advective  and  diffusive.  The  sandy  creek  channels  where 
the  majority  of  the  groundwater  (and  hence  nitrogen)  enters  the  marsh  make  up  rou^y 
13%  of  the  total  marsh  area,  while  the  vegetated  areas  where  it  is  utilized  cover  over  60% 
of  the  total  marsh  area.  The  disparity  in  these  two  percentages  illustrates  the  importance 
that  diffusion  must  have  in  determining  the  overall  nitrogen-utilization  capacity  of  the 
marsh  as  a  whole.  Consideration  of  only  advective  transport  may  be  sufficient  for 
determining  nitrogen  flux  to  the  coastal  ocean,  but  diffusive  transport  must  definitely  be 
considered  when  characterizing  processes  occurring  within  the  marsh  itself  Providing 
reliable  estimates  for  diffusivities  based  on  vegetation  densities  and  ambient  flows  may 
then  enable  us  to  fine  tune  the  marsh  to  improve  performance,  perhaps  by  altering 
vegetation  density  or  distribution  to  maximize  nutrient  removal  capability. 

These  diffusivities  depend  on  the  turbulence  levels  throughout  the  marsh,  with  higher 
turbulence  leading  to  greater  diffiisivity.  The  stem  size  and  population  density  of  marsh 
vegetation  can  be  a  key  determiner  of  turbulence  levels  in  the  marsh.  Although 
vegetation  coverage  is  generally  less  than  10%  in  natural  marshes  (Kadlec,  1990),  even 
this  minor  coverage  can  have  significant  impact  on  turbulence  (and  hence  diffiisivity) 
levels.  In  the  presence  of  adequate  flow,  addition  of  a  small  number  of  stems  increases 
turbulence  due  to  vortex  shedding.  This  effect  is  enough  to  counterbalance  the  decrease 
in  bed  production  associated  with  the  decrease  in  mean  flow  due  to  increased  drag  firom 


the  steins.  Increasing  the  stem  density  does  lead  to  a  reduction  in  turbulence  levels, 
however,  as  sufficient  flow  restriction  eventually  exists  to  limit  flow  enou^  to  prohibit 
vortex  shedding  (Nepf,  et  al,  1997,  Burke  and  Stolzenbach,  1983).  Based  on  a  numerical 
model  this  critical  stem  density  is  expected  to  be  around  1%  coverage,  with  turbulence 
levels  falling  off  at  both  higher  and  lower  stem  densities  (Burke  and  Stolzenbach,  1983). 
This  field  study  is  designed  to  explore  these  findings. 

It  should  be  noted  that  the  hydrodynamically-optimized  stem  density  discussed  above 
may  or  may  not  be  the  overall  ideal  stem  spacing  for  the  marsh  in  general.  For  example, 
it  has  been  noted  that  a  prolonged  increase  in  nutrient  loading  causes  Spartina 
alterniflora  plants  to  shift  from  a  short  to  a  tall  form  (V aliela  and  Teal,  1 979).  This 
change  in  plant  height  leads  to  a  concurrent  increase  in  average  plant  spacing,  potentially 
altering  the  size  and  types  of  species  able  to  use  the  marsh  for  protection,  foraging,  and 
other  purposes.  No  attempt  will  be  made  to  quantify  these  effects,  but  their  consideration 
may  be  critical  when  determining  what  the  truly  optimum  vegetation  density  is  for  a 
given  marsh. 

1.4  Other  Considerations 

As  discussed  above,  diffusivity  levels  observed  within  marsh  vegetation  are  dependent  on 
the  turbulence  structure  present.  While  the  goal  of  this  research  is  to  examine  diffusivity 
directly,  the  extension  to  turbulence  levels  may  also  be  useful  in  a  number  of  other  areas. 
For  example,  reduced  turbulence  and  flow  within  stands  of  vegetation  may  increase 
deposition  or  decrease  erosion  within  these  same  areas  -  a  fact  which  may  be  important  in 
the  study  of  sediment  transport  or  particulate  pollutant  deposition.  Reduced  turbulence 
may  also  be  an  important  determinant  of  food  supply  to  certain  marine  organisms,  since 
turbulent  mixing  may  be  critical  to  replenishing  the  food  supply  in  heavy  consumption 
areas  (Frechette,  et  al,  1989,  Butman,  et  al,  1994).  Additionally,  Anderson  and  Charters 
(1982)  observed  that  the  coastal  aquatic  plant  Gelidium  nudifrons  generates  turbulence  at 


flows  exceeding  6  cm/s  and  suggested  that  stem  size  is  adaptively  selected  to  optimize  the 
scale  of  turbulence  for  maximum  nutrient  uptake.  Other  researchers  have  related 
sediment-water  column  and  plant-water  column  exchange  of  oxygen,  phosphorus,  nitrate, 
and  ammonium  to  the  level  of  water  column  turbulence  (Koch,  1993,  Ziegler,  1993, 
Heuttel  and  Cost,  1992).  Since  turbulence  is  the  common  factor  in  all  these  concerns, 
relative  quantification  of  turbulence  levels  within  vegetation  stands  may  provide  useful 
insight  in  a  number  of  areas. 


2.  Study  Site  Description 


This  study  was  carried  out  in  Great  Sippewissett  Marsh,  a  tidal  salt  marsh  located  on 
Bu2zards  Bay  on  the  western  shore  of  Cape  Cod,  Massachusetts.  Great  Sippewissett 
Marsh  is  generally  characteristic  of  tidal  marshes  foimd  in  this  region,  specific 
characteristics  are  as  described  below. 

2.1  General  description 

Great  Sippewissett  Marsh  itself  occupies  483,800  m^,  with  approximately  44%  of  this 
area  being  classified  as  Tow  marsh’  and  20%  as  ‘high  marsh’  (Valiela  and  Teal,  1979). 
Short  and  tall  forms  of  Spartina  alterniflora  (cordgrass)  are  the  dominant  vegetation  in 
the  low  marsh  region,  while  Spartina  patens  and  Distichlis  spicata  dominate  the  high 
marsh.  The  remainder  of  the  marsh  is  broken  down  into  muddy  creek  bottoms,  sandy 
creek  bottoms,  and  algal  mats  -  with  each  being  present  in  roughly  equal  quantities. 

Seawater  fi:om  Buzzards  Bay  enters  and  leaves  the  marsh  through  a  single  entrance 
known  as  Sippewissett  Creek.  The  marsh  is  flooded  twice  daily  by  this  seawater,  with  the 
average  tidal  excursion  being  1 .6  meters  (Valiela  and  Teal,  1979).  Freshwater  enters  the 
marsh  through  a  combination  of  stream  and  groimdwater  inputs,  but  totals  only  about 
three  percent  of  the  saltwater  tidal  exchange  (Valiela  and  Teal,  1979).  Circulation  within 
the  marsh  is  therefore  almost  totally  controlled  by  tidally  driven  currents. 

Tidally  averaged  depths  in  Great  Sippewissett  Marsh  are  on  the  order  of  a  meter,  althou^ 
large  areas  of  the  marsh  are  very  shallow.  Figure  2-1  shows  an  overhead  view  of  Great 
Sippewissett  Marsh  with  the  study  area  and  other  prominent  features  called  out. 


2.2  Specific  study  conditions 


This  dye  release  experiment  was  carried  out  in  an  isolated  patch  of  Spartina  alterniflora 
located  in  the  Tow  marsh’  region  of  the  greater  marsh.  Figme  2-2  shows  this  patch  and 
the  area  immediately  surrounding  it.  The  patch  itself  was  roughly  3  meters  by  3  meters  in 
size  with  a  fairly  even  distribution  of  vegetation  over  this  area.  The  larger  area  around 
the  patch  had  a  relatively  smooth  and  flat  bed  {bed  being  the  lower  water  column 
boimdary)  and  was  exposed  during  low  tide  conditions,  allowing  both  for  ease  of 
equipment  setup  and  relatively  undisturbed  flow  conditions  during  measurements.  The 
study  area  was  located  immediately  to  the  south  of  a  major  west-to-east  creek  channel  in 
the  marsh  as  well,  so  the  currents  found  there  during  tidal  changes  are  believed  to  be 
among  the  highest  of  those  found  throu^out  the  marsh  (excluding  those  found  in  the 
creek  mouth  and  inlet  channel).  This  combination  of  strong  currents,  smooth  bed  and 
small  flow  resistance  seemed  to  make  this  an  ideal  site  for  the  dye  release  study.  Also, 
the  limited  extent  of  the  patch  encourages  more  uniform  flow,  i.e.  no  dead  end  effects 
caused  by  lack  of  an  outlet  for  the  water  flowing  into  the  vegetation.  Tidal  exclusions 
during  the  course  of  the  study  ranged  from  roughly  1.3  meters  early  on  to  0.9  meters  at 
the  conclusion. 


Figure  2-1;  Overhead  View  of  Great  Sippewissett  Marsh 


20 


21 


3.  Experimental  Techniques 


To  determine  diffusivities,  Rhodamine  WT  fluorescent  dye  was  continuously  injected  at 
the  midpoint  of  the  vegetation  patch,  and  the  resulting  fluorescence  profiles  were 
measured  using  a  flow-through  fluorometer  as  the  dye  plume  emerged  from  the  patch. 
These  profiles  were  then  used  to  estimate  diffusivities  for  the  vegetation  and  flow 
conditions  present.  This  chapter  discusses  the  equipment  used  to  accomplish  this 
injection  and  sampling,  operational  parameters  used  with  the  equipment,  specific 
sampling  techniques  used,  and  the  assumptions  used  in  gathering  and  analyzing  this  data. 

3.1  Equipment  Used 

Equipment  used  included  that  required  to  inject  die  Rhodamine,  to  collect  water  samples, 
and  then  to  analyze  these  samples.  Figure  3-1  shows  a  photograph  of  the  equipment  in 
the  field,  while  Figure  3-2  shows  a  schematic  drawing  of  the  general  layout.  More  details 
are  provided  in  the  paragraphs  which  follow. 

3.1.1  Fluorescent  Tracer  Addition 

The  tracer  solution  was  prepared  each  day  by  mixing  2-4  ml  of  a  20%  solution  of  the  dye 
in  10  liters  of  marsh  water.  The  resulting  mixture  was  then  injected  into  the  ambient  flow 
using  a  variable-stroke  positive  displacement  pump,  with  the  flow  rate  being  adjusted  as 
necessary  to  provide  adequate  tracer  concentrations  at  the  downstream  sampling  point. 
The  injection  rate  varied  from  day  to  day  depending  on  tidal  height  and  currents,  but 
generally  was  on  the  order  of  60  ml/minute.  This  allowed  a  sampling  window  of  up  to 
about  3  hours  without  depleting  the  dye  source. 

The  dye  injection  tube  itself  was  a  Vi”  ID  tygon  tube  secured  to  a  3/8”  diameter  threaded 
rod,  with  the  combination  (along  with  the  dye  bottle  and  the  injection  pump)  being 
mounted  on  a  0.5  meter  by  1  meter  table  constructed  from  %”  PVC  piping.  The  threaded 


Figure  3-2:  Schematic  Representation  of  Physical  Layout 
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rod  allowed  the  injection  height  to  be  easily  and  quickly  adjusted  vertically,  and  the 
dimensions  of  the  tube/rod  combination  were  selected  to  minimize  flow  disruption.  The 
mean  stem  diameter  of  vegetation  in  the  patch  was  foxmd  to  be  0.27”  (0.69  cm),  so  the 
tube/rod  combination  would  be  seen  much  like  another  stem  by  the  flow.  An  “L”  fitting 
on  the  free  end  of  the  tygon  tubing  directed  the  dye  discharge  parallel  to  the  ambient  flow, 
minimizing  any  mixing  due  to  shear  effects.  The  table  itself  was  oriented  with  its  long 
axis  perpendicular  to  the  ambient  current,  again  to  minimize  effects  from  the  table  legs  on 
the  turbulence  structure  near  the  dye  injection  point. 

3.1.2  Sample  Collection 

Water  was  collected  from  a  plane  located  roughly  1 .5  meters  downstream  from  the 
injection  point  by  means  of  a  sampling  pump  and  a  second  tygon  tube.  This  sampling 
tube  was  suspended  from  PVC  pipes  located  above  the  water’s  surface,  so  the  only 
intrusion  into  the  sampling  region  was  from  the  ‘A”  diameter  pipe  connected  to  the 
sampling  tube  and  used  to  adjust  it  vertically.  The  position  of  this  suspending  pipe  was 
adjustable  both  horizontally  and  vertically,  allowing  samples  to  be  gathered  throughout  a 
vertical  plane  immediately  downstream  from  the  vegetation  patch.  The  suspending  pipe 
itself  was  only  slightly  larger  diameter  than  a  standard  plant  stem  (again  minimizing  flow 
disruption),  and  the  sample  tube  itself  was  made  to  extend  roughly  5  cm  upstream  of  this 
suspending  pipe  to  further  minimize  any  disruptive  effects  near  the  sampling  point. 

The  sample  tube  was  ‘A”  ID  tygon  tubing,  and  the  sample  pump  drew  samples  at  a 
nominal  1 .5  gpm.  This  provided  for  a  velocity  at  the  sample  tube  inlet  of  roughly  75 
cm/second,  approximately  an  order  of  magnitude  higher  than  ambient  flows.  This 
condition  was  thought  to  be  unacceptable,  so  flie  end  of  the  sample  tube  was  fitted  with  a 
‘T’  fitting,  with  the  arms  of  the  “T”  trimmed  off  sli^tly,  to  minimize  this  problem.  This 
had  the  effect  of  roughly  halving  the  inlet  velocity  and  also  of  directing  the  suction  from 


the  sample  pump  perpendicular  to  ambient  current  flow.  This  condition  was  thou^t  to 
yield  more  acceptable  conditions  for  several  reasons: 

1 .  Distances  of  only  several  centimeters  from  either  end  of  the  “T”  are  required 
to  get  suction  flow  velocities  essentially  equal  to  ambient  values. 

2.  Dye  concentration  profiles  over  distances  such  as  these  (between  the  two 
inlets  plus  far  enough  into  the  fluid  for  suction  velocity  to  be  rou^y  ambient 
velocity)  should  be  well  approximated  linearly,  so  drawing  samples  from  both 
sides  may  tend  to  average  the  dye  concentration  very  close  to  the  true  value. 

While  this  condition  was  still  not  ideal,  it  was  deemed  acceptable  due  to  the  need  to 
sample  large  volumes  of  water  very  quickly  and  because  of  the  reasons  above. 

3.1.3  Sample  Analysis 

Samples  collected  with  the  above  apparatus  were  analyzed  in  situ  using  a  Turner  Designs 
Model  AU-10  fluorometer  in  continuous  flow  mode.  Water  was  drawn  through  the 
fluorometer  by  connecting  the  sample  tube  to  the  inlet  of  the  fluorometer's  sample 
chamber  and  pump  suction  to  the  outlet  of  the  sample  chamber.  A  background 
fluorescence  reading  was  obtained  at  the  beginning  of  each  day’s  sampling;  and  time, 
water  temperature  and  fluorescence  above  background  were  recorded  throughout  the 
sampling  period  using  the  internal  data  logger  installed  on  the  AU-10.  Further  details  on 
instrument  setup  and  data  collection  are  presented  in  the  following  sections. 

3.2  Data  Collection 

Data  collected  included  fluorescence  at  a  number  of  locations  in  the  vertical  plane  at  the 
downstream  boundary  of  the  vegetation  patch,  as  well  as  vegetation  density  and  ambient 
current  values.  Fluorescence  was  in  general  measured  in  one  vertical  profile  (same 


transverse  coordinate)  and  one  transverse  profile  (same  vertical  coordinate),  with 
additional  readings  obtained  as  possible.  More  detailed  information  on  this  data 
collection  is  presented  below. 

3.2.1  Sample  Duration 

Length  of  sampling  at  any  particular  coordinate  was  a  compromise  between  obtaining 
adequate  data  to  provide  a  representative  fluorescence  reading  at  that  point  and  obtaining 
readings  at  a  sufficient  number  of  these  coordinates  prior  to  significant  changes  in  the 
tidal  flow.  The  method  settled  upon  as  being  the  most  productive  involved  sampling  for 
one  minute  per  coordinate  for  one  full  vertical  and  one  full  transverse  profile,  and  then 
returning  to  take  longer  readings  at  each  coordinate  if  conditions  allowed.  In  general,  the 
one  minute  per  coordinate  snapshot  profiles  resulted  in  fluorescence  measurements  which 
were  more  precise  and  more  normally  distributed  than  the  longer  reading  profiles  due  to 
the  rapid  changes  in  tidal  flow. 

3.2.2  Sample  Location 

In  general,  a  vertical  profile  was  taken  at  what  appeared  to  be  the  horizontal  center  of  the 
plume,  followed  by  a  transverse  profile  at  the  vertical  coordinate  which  seemed  to  have 
the  mavimiim  averaged  fluorescence  reading.  These  averages  were  done  based  on  the 
visual  display  of  the  AU-10,  and  so  were  not  always  precise  as  shown  by  later,  more 
detailed,  data  analysis.  Readings  for  any  given  profile  were  normally  taken  at  10  cm 
increments,  althou^  finer  and  more  coarse  profiles  were  also  taken  as  time  and 
conditions  allowed. 

3.2.3  Vegetation  Density 

Vegetation  density  was  determined  by  randomly  casting  a  0. 1  m  firame  onto  the 
vegetation  island  and  coimting  the  number  of  stems  within  the  firame.  This  measurement 
was  repeated  five  times  for  each  vegetation  density  to  be  measured,  with  the  density  used 


for  further  calculations  being  the  mean  of  the  five  individual  measurements.  These 
measurements  were  fibrst  taken  at  the  natural  density,  then  stems  were  selectively  (and 
randomly)  removed  in  a  series  of  steps  to  gradually  reduce  vegetation  density  to  zero.  In 
each  step,  roughly  one-half  of  the  remaining  stems  were  removed.  This  removal  was 
done  by  trimming  off  the  stems  3-4  centimeters  above  the  bottom,  allowing  the  stem  to 
be  removed  without  killing  the  plant.  (Full  vegetation  density  has  already  been  restored 
due  to  stem  regrowth.)  A  representative  sample  of  roughly  100  stems  from  the  first 
cutting  was  measured  for  diameter  at  the  base  of  the  stem;  mean  stem  diameter  was  found 
to  be  0.69  ±  0.03  cm.  This  value  was  then  used  to  compute  percent  vegetation  coverage, 
defined  as  the  percentage  of  the  bed  area  actually  occupied  by  plant  stems  and  designated 
as  P.  For  example,  a  stem  density  of  200  stems/m^  gives  a  total  stem  area  of  74.8  cm^ 
(200*;t*(0.69/2)^),  resulting  in  a  percent  vegetation  coverage  of  0.75%.  This  technique 
does  not  consider  the  effects  of  any  leaves  attached  to  the  stem  or  any  variation  in  stem 
size  with  depth,  but  since  these  effects  should  (at  least  to  the  first  order)  be  proportional 
to  stem  density,  it  provides  a  reasonable  measure  of  vegetation  coverage.  Table  3-1 
shows  the  vegetation  densities  and  percent  stem  coverages  found  in  this  study. 


Table  3-1:  Vegetation  Densities  and  Percent  Coverages 


Vegetation 

Density 

Stems/0.1  ni2 

(each  of  5  trials) 

Stems/m2 

(mean) 

Percent 

Stem 

Coverage 

High 

40, 36,  37, 37, 35 

370 

1.38 

Medium 

19, 17,21,21,20 

196 

0.73 

Low 

9,8,11,10,10 

96 

0.36 

Zero 

stubble  only 

0 

0 
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3.2.4  Determination  of  Currents 


Currents  were  measured  by  injecting  a  parcel  of  dye  at  roughly  mid-depth  and  timing  its 
transit  across  a  known  distance.  This  measurement  was  repeated  for  a  number  of  trials  at 
several  different  times  (roughly  15  minute  intervals)  over  the  course  of  a  given  sample 
profile,  with  the  individual  readings  at  a  given  time  being  averaged  to  give  the  mean 
current  for  that  time.  The  mean  sample  time  for  all  coordinates  used  in  a  given  profile 
was  also  obtained,  and  a  current  for  this  time  was  found  via  linear  interpolation  between 
the  measured  ciurent  values  on  either  side  of  the  time  in  question.  This  cinrent  was  then 
used  in  all  further  calculations  for  that  particular  profile. 

3.3  Fluorometer  Setup 

The  Tmner  Designs  Model  AU-10  fluorometer  used  in  this  study  has  several  user- 
definable  parameters  for  customizing  the  data  gathered.  Settings  used,  and  the  rationale 
behind  these  choices,  are  discussed  in  the  sections  below. 

3.3.1  Data  Averaging  Interval 

The  visual  display  on  the  AU-10  is  updated  with  a  new  fluorescence  reading  every  one 
second;  the  value  is  recorded  in  the  internal  data  logger  is  either  these  instantaneous 
values  or  some  time-average  of  these  readings.  Readings  were  logged  every  one  second 
throughout  this  study  to  maximize  data  availability,  so  the  data  averaging  interval  was 
selected  as  one  second.  This  in  effect  records  the  instantaneous  values,  however,  since 
only  one  reading  is  averaged  to  obtain  the  recorded  fluorescence. 

3.3.2  Time  Constant 

The  fluorometer  also  has  an  internal  time  constant  setting  which  works  to  reduce 
variability  in  fluorescence  readings.  As  discussed  above,  the  fluorometer’ s  display  is 
updated  every  second,  but  the  AU-10  is  actually  sampling  internally  at  10  Hertz.  These 
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internal  readings  are  then  averaged  (using  the  time  constant  setting)  to  give  the  value 
displayed  on  the  visual  display  and  sent  to  the  internal  data  logger.  A  time  constant  of  1 
second  was  used  throughout  this  study.  This  setting  causes  the  AU-10  to  use  the  previous 
1 .6  seconds  worth  of  internal  readings  to  determine  the  value  to  be  displayed/sent  to  the 
data  logger,  and  is  die  shortest  time  constant  available  without  modifying  the  instrument. 
Alternate  time  constants  available  are  2, 4,  or  8  seconds;  these  use  the  previous  3.2, 6.4, 
and  12.8  seconds  worth  of  internal  readings,  respectively,  when  determining  the 
fluorescence  value  displayed  and  sent  to  the  logger.  Using  the  1  second  time  constant 
then  allowed  for  maximum  instrument  response  and  minimized  any  internal  smoothing  of 
the  data. 

3.3.3  Temperature  Compensation 

Fluorescence  decreases  with  increasing  temperature  in  accordance  with  Equation  3-1. 

F  =  Fo*  ^  Equation  3-1 

where  F  =  fluorescence  reading  at  a  given  temperature 

Fo  =  fluorescence  reading  at  the  chosen  reference  temperature 

n  =  the  temperature  exponent  for  the  dye  in  use 

to  =  the  chosen  reference  temperature 

o 

The  reference  temperature  was  chosen  as  25  C  for  this  study,  and  the  temperature 
exponent  was  set  at  -0.027  for  Rhodamine  WT  as  suggested  in  the  literature  (Smart  and 
Laidlaw,  1977).  Experiments  conducted  using  a  closed-loop  recirculation  system  showed 
that  this  compensation  scheme  removed  roughly  75%  of  the  effect  of  temperature  on 
fluorescence  readings,  but  left  a  variation  of  2%  of  the  fluorescence  reading  per  degree  C. 
This  is  not  as  exact  a  compensation  as  might  be  hoped  for,  but  since  temperature 
variations  over  a  sampling  period  were  at  most  2'-3°C,  the  maximum  effect  on  relative 


fluorescence  readings  should  be  rou^y  5%.  Temperature  changes  over  any  individual 
profile  were  never  more  than  0.5°C,  however,  so  errors  due  to  temperature  compensation 
effects  should  be  1%  over  any  profile.  This  error  is  well  within  the  uncertainty  associated 
with  any  mean  fluorescence  value,  and  so  may  considered  negligible.  One  must  be 
conscious  of  this  potential  disparity,  however,  when  comparing  absolute  fluorescence 
readings  at  differing  temperatures  rather  than  relative  readings  across  a  single  profile. 

3.4  Simplifying  Assumptions 

A  number  of  simplifying  assumptions  were  made  in  this  study  due  to  equipment  and 
personnel  availability.  Specific  assumptions  are  discussed  below. 

3.4.1  Uniform  current  profile 

Currents  within  the  patch  of  vegetation  were  measured  by  injecting  a  parcel  of  dye 
approximately  midway  between  the  water  surface  and  the  bed  and  timing  its  transit  across 
a  measured  distance.  This  injection  was  done  only  at  one  height  within  the  canopy,  so  the 
current  was  actually  measured  only  at  this  one  height.  Vertical  variations  in  plant 
morphology  (and  thus  drag)  might  be  expected  to  create  some  vertical  variation  in  the 
current  profile,  although  studies  have  shown  (Leonard  and  Luther,  1995,  Gambi,  et  al, 
1990,  Shi,  et  al,  1996,  El-Hakim  and  Salama,  1992)  this  variation  within  vegetation 
canopies  to  be  relatively  small  -  generally  on  the  order  of  10%  of  the  mean  velocity. 
Therefore,  the  simple  single-point  measurement  obtained  should  give  a  reasonable 
approximation  for  die  depth  averaged  velocity. 

The  current  was  also  assumed  to  be  imiform  throughout  the  spatial  extent  of  the 
vegetation  patch.  This  is  likely  to  be  true  except  very  near  the  edges  of  the  patch  (due  to 
boundary  effects),  and  since  no  injections  or  measurements  were  made  near  the  edges  of 
the  patch,  seemed  to  be  a  good  assumption.  This  uniform  three-dimensional  current 


assumption  then  removes  the  effect  of  shear  dispersion,  leaving  turbulent  diffusion  effects 
as  the  major  parameter  being  measured. 

3.4.2  Boundary  Effects  Negligible 

Since  the  goal  was  to  isolate  the  contributions  of  the  plant  stems  to  turbulence,  all  other 
effects  had  to  be  minimized.  These  included  the  influence  of  nonuniform  conditions 
present  near  the  edges  of  the  vegetation  patch,  turbulence  introduced  by  the  bed, 
turbulence  structure  carried  into  the  patch  with  the  current,  and  wind  and  wave  effects. 
The  first  factor  was  discussed  above  and  was  assumed  negligible  since  no  injections  or 
measurements  were  made  near  the  patch  boxmdaries.  All  measurements  were  in  fact 
taken  at  least  0.5  meters  firom  the  edge  of  the  patch,  the  best  that  could  be  done  with  a 
patch  of  this  size. 

To  minimize  contribution  firom  the  bed,  dye  injections  were  made  approximately  midway 
between  the  surface  and  the  bed,  and  measurements  were  taken  longitudinally  close 
enough  to  the  injection  point  so  that  the  plume  had  not  yet  grown  sufficiently  in  the 
vertical  to  feel  the  presence  of  the  bed.  It  is  probable  that  bed-generated  turbiilence  was 
present  at  the  height  of  the  plume,  but  stem-generated  turbulence  generally  dominates  the 
contribution  firom  the  bed  (Nepf,  et  al,  1997),  at  least  at  higher  stem  densities.  At  lower 
stem  densities  this  assumption  may  not  be  as  valid,  since  reduced  stem-generated 
turbulence  may  cause  bed-generated  turbulence  to  become  more  significant.  These 
findings,  coupled  with  the  mid-water  column  injection  and  close-in  sampling,  make  it 
reasonable  to  neglect  bed-generated  turbulence  within  the  vegetation  patch,  at  least  as  a 
first  approximation. 

Turbulence  structure  carried  into  the  vegetation  patch  was  also  a  concern.  The  site  was 
selected  to  minimize  this  contribution,  with  approximately  25  meters  of  relatively  straight 
flow  channel  with  a  smooth  and  flat  bed  existing  both  up-  and  downstream  of  the 
vegetation  patch.  This  minimized  any  flow  disruptions  and  allowed  for  data  gathering  on 
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both  incoming  and  outgoing  tides.  Additionally,  since  the  dye  injection  point  was 
essentially  in  the  center  of  the  vegetation  patch,  approximately  1 .5  meters  of  vegetation 
provided  a  zone  of  adjustment  before  the  flow  reached  the  injection  point.  These 
characteristics  combined  to  minimize  the  effects  of  externally  induced  turbulence  within 
the  vegetation  patch. 

Wind  and  wave  effects  were  not  specifically  accoimted  for,  althou^  fliere  also  were  no 
specific  days  which  seemed  significantly  more  (or  less)  windy  or  wavy  than  any  other 
days.  These  effects  are  then  assumed  to  be  relatively  insignificant,  althou^  no  hard 
evidence  was  analyzed  to  support  this  conclusion. 

3.4.3  Fluorescence  Changes  Negligible 

Changes  in  the  fluorescence  conditions  could  occur  during  the  course  of  die  concentration 
measurements.  These  changes  could  be  due  either  to  changes  in  the  water  temperature 
affecting  fluorescence  readings  or  to  changes  in  the  composition  of  the  marsh  water. 
Fluorescence  readings  change  as  temperature  changes.  Temperature  variations  over  the 
full  tidal  cycle  were  routinely  seen  to  be  roughly  3°C,  but  the  range  over  any  related  series 
of  measurements  (taken  over  a  period  of  at  most  30  minutes)  would  then  be  less  than  half 
a  degree.  This,  coupled  with  the  temperature  compensation  feature  of  the  fluorometer 
used,  made  this  effect  negligible.  Fluorescence  readings  could  also  change  with  water 
composition,  due  to  either  background  fluorescence  changes  or  changes  in  salinity. 
Measurements  of  background  fluorescence  over  a  full  tidal  cycle  showed  negligible 
change.  These  results  were  then  extrapolated  to  later  days  and  this  effect  was  assumed 
not  to  be  important  throughout  the  study.  Changes  in  fluorescence  readings  due  to 
salinity  are  more  complex.  There  is  a  marked  decrease  in  measured  fluorescence  from 
Rhodamine  WT  as  salinity  increases,  with  the  most  significant  breaking  point  being 
around  0.03  M  NaCl  (2  ppt  salinity)(Smart  and  Laidlaw,  1977).  The  decrease  is  also 
quite  time  dependent,  however,  with  other  studies  suggesting  minimal  initial  effect 


followed  by  major  changes  as  the  sample  solution  stands  for  several  days.  For  example, 
0.6  M  NaCl  (35  ppt  salinity)  causes  a  50%  reduction  in  fluorescence  over  a  27  hour 
period  (Smart  and  Laidlaw,  1977).  A  series  of  comparative  measurements  between 
standards  prepared  xxsing  distilled  water  and  several  marsh  water  samples  showed  a 
decrease  in  fluorescence  of  4-8%  for  the  marsh  water  samples  over  the  distilled  water 
samples.  This  demonstrated  minor  effect  over  a  range  of  possible  salinities,  coupled  with 
the  relatively  short  elapsed  time  during  any  series  of  related  measurements,  led  to  this 
effect  being  assumed  as  negligible  in  this  study. 

3.4.4  Conservative  Tracer 

The  Rhodamine  WT  dye  used  as  the  tracer  in  this  experiment  was  considered  to  behave 
conservatively.  This  assumption  took  into  account  the  dependence  on  salinity  as 
discussed  above  as  well  as  adsorptive  and  photochemical  losses.  Rhodamine  WT  is  one 
of  the  most  conservative  tracers  available,  with  minimal  adsorptive  or  decay  losses, 
particularly  at  low  concentrations.  Several  studies  have  measured  overall  losses  at  less 
than  5%  for  this  tracer  (Smart  and  Laidlaw,  1977),  so  these  losses  were  not  considered  in 
this  study.  During  the  course  of  the  study,  no  evidence  of  adsorption  to  the  vegetation 
was  found,  although  there  did  appear  to  be  some  adsorption  to  sediments  well  outside  the 
study  area.  This  was  based  on  a  distinct  pink  color  present  in  the  sediment  for  several 
days  after  the  dye  release. 


4.  Analysis  Techniques 


This  chapter  presents  the  techniques  used  in  analyzing  the  data  gathered  as  described  in 
the  previous  chapter.  Example  analyses  and  calculations  are  provided  using  a 
representative  profile. 


4.1  Selection  of  Data 

The  sampling  port  was  moved  between  locations  at  set  times,  but  since  this  movement 
was  not  instantaneous  and  since  the  sample  pump  was  continuously  drawing  samples, 
some  of  the  data  points  gathered  were  not  truly  representative  of  the  conditions  at  the 
intended  sample  point.  Points  collected  during  these  transitions  between  sampling  points 
were  cut  from  the  data  record  based  on  the  following  discussion. 


The  sample  line  consisted  of  roughly  20  feet  of  Vi  inch  ID  tygon  tubing  with  a  small  filter 
located  about  halfway  between  the  sample  inlet  and  the  fluorometer.  This  volume,  along 
with  the  small  volume  inside  the  fluorometer’ s  sample  chamber,  totaled  approximately 
one  liter  and  represented  the  volume  to  be  purged  between  samples.  Since  the  sample 
pump  had  a  nominal  capacity  of  1.5  gpm  (5.7  liter/minute),  this  represented  rnngbly  10 
seconds  of  pranping.  The  first  thirty  seconds  of  data  fi-om  each  sample  point  were 
actually  discarded,  however,  to  allow  for  file  time  delay  in  shifting  the  sample  line  and  to 
allow  for  adequate  purging  and  instrument  stabilization.  This  allowed  for  purging 
roughly  three  line  volumes  before  gathering  data  which  were  actually  used  in  the  analysis; 
these  discarded  points  were  recorded  and  are  available  for  use  if  needed,  however. 


Figure  4-1  shows  the  full  horizontal  data  series  for  the  representative  profile,  while  Figure 
4-2  shows  the  sections  of  this  profile  used  to  calculate  fluorescence  at  each  coordinate. 
The  coordinates  listed  in  Figure  4-2  are  transverse  and  vertical  displacement  fi-om  an 
arbitrary  reference  point.  As  discussed  earlier,  the  first  30  data  points  at  each  coordinate 
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Profile  Point  3  -  (120-40) 


Profile  Point  4 -(130^0) 


Figure  4-2:  Data  Used  for  Individual  Profile  Points 
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were  not  analyzed,  and  all  remaining  data  points  were  used.  Each  coordinate  in  this 
particular  profile  then  had  30  data  points  which  were  used  in  analysis. 

4.2  Analysis  Methods 

Remaining  data  points  for  each  sample  location  were  then  averaged  to  get  the  mean 
fluorescence  reading  at  that  location.  The  standard  deviation  about  this  mean  was  also 
calculated,  and  the  error  m  the  calculated  mean  value  was  estimated  using  Equation  4-1 . 


§  =  Equation  4-1 

where  dp  is  the  error  in  the  mean  fluorescence  and  n  is  the  number  of  data  points  used  in 
finHing  the  mean.  This  technique  assumes  that  the  individual  data  points  used  are 
independent,  random  samples  -  this  is  not  absolutely  true  here,  but  this  method  still  gives 
reasonable  estimates  on  the  uncertainty  associated  with  the  mean  as  n  becomes  large. 
Table  4-1  shows  the  mean  fluorescence  and  associated  uncertainty  for  each  of  the  points 
in  the  representative  profile;  the  standard  deviation  is  also  tabulated  for  use  later. 


Table  4-1:  Mean  Fluorescence  Values  for  Profile  Points 


Data  Point 

Coordinates 

Mean 

Fluorescence 

(PPb) 

Uncertainty 
in  Mean 
Fluorescence 
(ppb) 

Standard 

Deviation 

for 

Fluorescence 

1 

100-40 

37.55 

4.90 

13.7 

2 

110-40 

296.83 

2.56 

7.15 

3 

120-40 

55.20 

2.44 

6.81 

4 

130-40 

0.60 

0.19 

0.52 

Once  the  average  fluorescence  values  were  determined,  the  distribution  of  mean 
concentrations  in  space  was  evaluated  using  two  different  methods  as  discussed  below. 
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4.2.1  Method  of  Moments 


The  first  method  followed  the  Method  of  Moments  analysis  presented  in  Fischer,  et  al., 
(1979,  p.  39-40),  and  used  Equations  4-2  throu^  4-6.  A  linear  interpolation  was  used  to 
approximate  the  fluorescence  profile  between  adjacent  data  points,  and  a  linear 
extrapolation  to  zero  fluorescence  was  done  using  the  first  two  and  the  last  two  data 
points  in  the  profile. 


oo 


M0=  ^C{x,t)dx 

— oo 

Equation  4-2 

oo 

M\=  ^xC{x,t)dx 

— oo 

Equation  4-3 

00 

M2=  ^x^C{x,t)dx 

—00 

Equation  4-4 

x  =  Ml/MO 

Equation  4-5 

00 

cr"  =  \{x'fC{x\t)dx'IMO  =  (M2  /  MO) 

Equation  4-6 

—00 


The  value  of  x  in  Equations  4-3  and  4-4  above  is  really  distance  fi:om  the  plume  center  of 
mass,  and  since  coordinates  were  measured  relative  to  a  fixed  point  some  conversion  had 
to  occur.  Equations  4-2, 4-3,  and  4-5  were  first  used  to  determine  the  coordinate  of  the 
center  of  mass  of  the  profile,  then  this  value  was  subtracted  from  all  x  coordinates  to  give 
new  coordinates  (x*)  relative  to  the  center  of  mass.  This  technique  set  x  (and  so  Ml)  to  0, 
and  x'  values  ranged  over  both  positive  and  negative  values.  Equations  4-2  and  4-4  were 
then  properly  evaluated  using  these  new  x'  coordinates,  and  finally  Equation  4-6  was  used 
to  find  the  variance  of  the  profile  using  the  values  found  for  M2  and  MO. 


Figure  4-3  shows  the  four  profile  points  used  in  the  representative  analysis  as  well  as  the 
two  extra  points  found  by  extrapolating  fluorescence  to  zero  (the  extra  point  at  the  right 
end  lies  at  130.1,  and  is  basically  the  same  as  the  last  real  point);  the  dashed  line 
connecting  the  points  shows  the  linear  approximation  to  the  normal  profile.  Analysis  as 
described  above  gave  a  plume  center  of  mass  at  1 1 1 .1  cm,  a  standard  deviation  of  5.73 
cm,  and  a  variance  of  32.80  cm  . 


350 

S'  300 
^  250 
I  200 
I  150 
§  100 
^  50 

0 

80  90  100  no  120  130  140 

Horizontal  Coordinate  (cm) 


Figure  4-3:  Method  of  Moments  Analysis  -  Representative  Profile 
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4.2.2  Exponential  Fit 


The  fluorescence  profiles  were  expected  to  be  normally  distributed,  so  the  second  form  of 
analysis  involved  fitting  a  standard  normal  equation  to  the  measured  data  points.  This 
standard  normal  equation  is  given  by  Equation  4-7,  the  solution  was  obtained  by 
linearizing  it  as  in  Equation  4-8  and  then  applying  a  weighted  linear  least  squares 
technique.  Uncertainties  in  fluorescence  readings  were  converted  to  uncertainties  on  a 
logarithmic  basis  using  Equation  4-9,  then  weights  for  the  individxial  points  were 
computed  using  Eq\iation  4-10.  (Taylor  (1982,  p.  150)  describes  the  weighting  of 
individual  points  by  the  inverse  of  the  square  of  their  standard  deviation.)  This  technique 
resulted  in  those  points  with  greater  relative  uncertainty  in  fluorescence  receiving  less 
weight  than  more  certain  points  when  doing  the  least  squares  analysis.  The  x  value 
found  using  the  Method  of  Moments  was  used  as  a  starting  point  for  this  exponential  fit 
analysis,  with  values  ±5  cm  from  this  starting  point  being  checked  m  0.1  cm  increments 
to  examine  quality  of  fit.  The  3c  value  over  this  range  giving  the  best  linear  correlation 
coefficient  for  the  least  squares  line  was  chosen  as  the  best  center  of  mass,  and  the 
constants  A  and  B  determined  for  this  point  were  used  in  further  analyses. 

F  =  Ae  Equation  4-7 

ln(F)  =  -B(x  -  xy  +  ln(.4)  Equation  4-8 

=<7pl  F  Equation  4-9 

)  Equation  4-10 

Analysis  of  the  representative  profile  using  this  second  technique  resulted  in  a  center  of 
mass  at  1 10.1  cm,  with  this  value  giving  an  value  of  0.991  for  the  least  squares  fit. 
Weightings  used  for  points  1  through  4  were  7.5, 1724.3, 65.6,  and  1.3,  and  best  values 
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for  ^  and  B  were  5.69  and  0.017,  respectively.  Figure  4-4  shows  the  individual  profile 
points  and  the  best  fit  exponential  curve  found  using  this  method. 


Figure  4-4:  Best  Fit  Exponential  Curve  -  Representative  Profile 


4.3  Calculation  of  Diffusivity 

Diffiisivity  for  a  plume  at  steady  state  can  be  found  using  Equation  4-11,  with  u  being  the 
mean  current  velocity  and  /  being  the  longitudinal  distance  fi:om  the  source.  The  variance 
)  is  the  direct  ouq)ut  of  the  Method  of  Moments  analysis,  so  calculation  of  diffusivity 
is  straightforward  in  this  case.  The  value  of  B  foimd  in  the  exponential  fit  method  is 
equivalent  to  l/2a  ,  so  this  calculation  is  also  straightforward  after  making  the  conversion 
fi'om  B  to  variance. 


Equation  4-11 
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The  longitudinal  displacement  Q)  for  the  representative  profile  was  127  cm.  As  discussed 
in  section  3.2,  currents  for  profile  mean  times  were  determined  by  linearly  interpolating 
between  currents  measured  at  roughly  15  minute  intervals  throughout  each  day’s  study. 
Equation  4-12  describes  the  interpolation  between  points  ti  and  t2,  while  Figure  4-5 
shows  the  measured  currents  for  the  day  of  the  representative  profile.  Data  for  profile 
points  1-4  were  taken  between  2:25  and  2:29,  giving  a  mean  time  of  2:27  for  this  profile; 
interpolation  between  the  first  two  points  in  the  current  profile  gives  a  current  of  2.07 
cm/s  for  that  time.  This  current  was  then  used  in  evaluating  diffiisivities. 


u{t)  =  u{t^)  + 


(/-?])  =  «(/, )  +  {slop^{t  -  r, ) 


Equation  4-12 


Values  obtained  for  diffusivity  were  either  transverse  (horizontal)  or  vertical  depending 
on  the  profile  being  analyzed.  Transverse  values  were  designated  as  Dx,  while  vertical 
values  were  designated  as  Dy.  For  this  profile,  Dx  was  foimd  to  be  0.27  cm^/s  using  the 
Method  of  Moments  and  0.24  cm  /s  using  the  exponential  fit  technique. 
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Figure  4-5:  Representative  Current  Profile 
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4.4  Evaluation  of  Uncertainty 


Uncertainties  associated  with  final  and  intermediate  values  were  determined  using  the 
technique  described  by  Taylor  (1982).  Equation  4-13  gives  this  relation  for  the 
propagation  of  uncertainty  in  a  function  of  several  variables,  this  method  was  used  in 
evaluating  uncertainties  in  currents,  plume  variance,  and  dififiisivity  as  discussed  below. 


Equation  4-13 


4.4.1  Uncertainty  in  Current 

As  discussed  earlier,  currents  at  any  given  time  were  obtained  by  linearly  interpolating 
between  two  measured  values  as  shown  by  Equation  4-12.  The  uncertainty  in  this 
interpolated  current  is  then  related  to  the  uncertainties  in  the  two  measured  currents 
according  to  Equations  4-14  and  4-15.  Since  times  were  recorded  to  the  nearest  minute, 
the  uncertainty  in  the  time  was  set  at  1  minute  for  this  analysis. 


For  the  representative  profile,  uncertainty  in  the  interpolated  current  (2.07  cm/s  at  2:27) 
was  found  to  be  +  0.19  cm/s,  based  on  interpolation  between  the  1.12  ±  0.09  cm/s  value 
at  2:1 1  and  the  2.55  ±  0.22  cm/s  value  at  2:35. 
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4.4.2  Uncertainty  in  Plume  Variance 


4.4.2. 1  Method  of  Moments 

Since  in  this  case  variance  is  estimated  by  trapezoidal  integration,  the  resulting 
uncertainty  is  due  both  to  uncertainty  in  the  individual  data  points  as  well  as  the  errors 
introduced  by  approximating  the  smooth  fluorescence  profile  as  a  series  of  trapezoids. 
This  second  effect  is  difficult  to  quantify  reliably,  since  its  effects  are  likely  to  vary  from 
profile  to  profile.  Additionally,  the  magnitude  of  this  contribution  is  also  likely  much  less 
than  that  from  uncertainty  m  the  data  itself,  so  efforts  focused  on  estimating  the 
uncertainty  due  to  the  data. 

To  do  this  analysis,  each  data  point  was  assigned  three  values  -  die  mean  fluorescence  and 
the  mean  plus  and  minus  the  uncertainty  in  fluorescence.  The  trapezoidal  integration  was 
then  done  using  every  possible  permutation  of  these  values  for  all  data  points  in  the 
profile.  The  mean  of  the  resulting  variance  values  was  then  found,  and  the  uncertainty  in 
this  mean  was  set  at  plus  and  minus  twice  the  standard  deviation  of  these  variances.  This 
technique  results  in  a  sort  of  worst  case  error  estimate,  but  seems  like  a  reasonable 
approach  given  the  sparsity  of  the  data  and  their  varying  imcertainties.  This  conservative 
approach  should  also  help  to  account  for  any  errors  introduced  by  the  trapezoidal 
integration  scheme  and  not  accounted  for  in  the  formal  error  analysis.  For  the 
representative  profile,  plume  variance  ranged  from  31.80  to  33.82  with  a  mean  of  32.80 
and  a  standard  deviation  of  0.51,  so  plume  variance  was  set  at  32.80  ±  1.02  cm  /s. 

4.4.2.2  Exponential  Fit 

In  this  case,  uncertainties  in  the  data  were  accounted  for  in  calculating  the  exponential 
constant  B.  The  variance  in  fi  is  given  by  Equation  4-16  (Taylor,  1982),  and  the 
uncertainty  in  B  was  set  at  two  times  the  standard  deviation  (giving  the  95%  confidence 
level)  found  from  this  equation.  This  uncertainty  in  B  was  converted  to  uncertainty  in 


profile  variance  by  applying  Equation  4-13  to  the  relation  between  B  and  cP"  to  obtain 
Equation  4-17.  For  the  representative  profile,  B  was  found  to  be  0.0170  ±  0.0031  cm'^, 
resulting  in  a  profile  variance  of 29.36  ±  5.29  cm^  using  this  method. 


Equation  4-16 


Equation  4-17 


4.4.3  Uncertainty  in  Diffusivity 

Equation  4-18  provides  the  uncertainty  in  the  difiusivity  based  on  the  imcertainties 
calculated  above,  and  was  obtained  by  applying  Equation  4-13  to  Equation  4-11. 
Uncertainty  in  /  was  set  at  2  cm  for  all  analyses,  uncertainty  in  all  other  parameters  was 
calculated  for  each  profile  individually.  Horizontal  diffusivity  was  found  to  be  0.27  ± 
0.03  cm^/s  using  the  Method  of  Moments  approach,  and  0.24  ±  0.05  cvP/s  using  the 
exponential  fit  technique  for  the  representative  profile. 


SD  = 

Analysis  for  all  other  profiles  is  similar  to  that  presented  for  the  representative  profile. 
Results  firom  these  analyses  are  summarized  in  the  next  chapter.  More  detailed 
information;  such  as  profile  point  coordinates,  fluorescence,  and  currents;  is  found  in 
Appendix  A  at  the  end  of  this  thesis. 
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Equation  4-18 


4.5  Reliability  of  Data 


The  only  remaining  factor  influencing  the  viability  of  the  results  obtained  is  the  reliability 
of  the  data  used  to  generate  them.  As  mentioned  earlier,  the  sampling  routine  for  this 
experiment  became  a  balance  between  gathering  sufficient  data  at  each  individual  point 
and  gathering  a  sufficient  number  of  these  points  during  a  given  tidal  excursion  to  give  a 
representative  profile.  The  tide  comes  in  and  out  of  Great  Sippewissett  Marsh  very 
quickly,  generally  allowing  only  rou^y  30  minutes  of  data  collection  on  each  tide.  The 
technique  that  seemed  to  produce  the  best  profiles  given  the  time  available  entailed  taking 
one  minute  of  data  at  each  point,  creating  a  snapshot  of  the  profile  over  a  3  -  6  minute 
period.  This  routine  did  by  no  means  encompass  the  variability  of  the  fluorescence 
readings  at  these  individual  points,  however,  and  even  these  snapshots  are  subject  to  the 
variabilities  in  conditions  over  the  time  needed  to  produce  them.  By  visual  observation,  it 
would  be  hard  to  say  if  the  plume  ever  did  really  reach  equilibrium,  since  rising  and 
falling  tides  continually  change  the  hydrodynamic  conditions  which  govern  the  flow 
through  the  patch.  So  while  this  sampling  routine  may  have  been  sufficient  to  resolve 
some  of  the  smaller  time  scale  diffiisivity,  longer  time  scale  events  will  not  be  seen  by 
this  technique.  As  discussed  earlier.  Ibis  may  have  been  the  cause  of  the  lower  than 
expected  diffiisivities  seen  for  the  vegetation-firee  profiles.  These  longer-term 
meanderings  and  alterations  of  the  ambient  flow  and  larger  scale  turbulence  events  could 
contribute  significantly  to  diffusion  in  the  tidal  marsh,  however,  so  proper  assessment  of 
their  contributions  would  be  necessary  to  more  acciirately  quantify  predicted  diffiisivities. 


5.  Results  and  Discussion 


5.1  General  Observations 

In  general,  measured  values  for  transverse  and  vertical  diffiisivities  were  essentially  equal 
to  each  other  and  behaved  similarly  across  all  ranges  of  vegetation  density  and  ambient 
flow.  The  magnitudes  of  these  values  were  also  on  the  order  of  those  previously 
observed  in  natural  settings  and  laboratory  studies  (Worcester,  1995,  Nepf,  et  al,  1997). 
Furthermore,  although  in  many  cases  there  is  significant  uncertainty  in  the  calculated 
dififusivities,  the  overall  trend  shows  a  roughly  linear  increase  in  diffusivity  with  ambient 
flow.  This  result  is  not  unexpected,  since  turbulence  intensity  (and  hence  turbulent 
diffusivity)  should  increase  as  flow  increases.  Diffusivities  were  not  found  to  be 
dependent  on  vegetation  density,  and  diffusivities  for  vegetated  conditions  compared  very 
well  with  those  for  non-vegetated  conditions.  More  specifics  on  these  observations  and 
some  thoughts  on  their  causes  are  presented  in  what  follows. 

Table  5-1  and  Table  5-2  provide  the  results  of  this  study  in  tabular  form.  The  column 
headings  for  these  tables  are  as  follows:  P  is  percent  vegetation  coverage,  del  x  refers  to 
the  uncertainty  in  x,  D-MOM  is  the  diffusivity  calculated  by  the  Method  of  Moments,  and 
D-EF  is  the  diffusivity  calculated  by  the  exponential  fit  method.  Each  run  is  also 
assigned  a  three  character  designator  identifying  it  as  either  a  horizontal  or  vertical 
profile,  a  sequential  ntimber,  and  specifying  either  low,  medium^  hi^,  or  zero  vegetation 
density.  For  example,  profile  H-4L  is  the  fourth  horizontal  profile  in  low  vegetation 
density.  Several  profiles  gave  anomalous  or  questionable  diffusivity  results,  these  are 
discussed  further  in  the  notes  which  follow  Tables  5-1  and  5-2. 


Table  5-1:  Transverse  Diffusivity  Results  for  Field  Study 


Run  # 

P 

BiWi 

delu 

(cm/s) 

Re-d 

D-MOM 

(cm''2/s) 

delD 

(cm''2/s) 

D-EF 

(cm^2/s) 

delD 

(cm''2/s) 

Depth 

(cm) 

H-IL 

0.36 

7.40 

0.43 

511 

0.97 

0.12 

0.82 

0.07 

42 

H-2L 

0.36 

4.62 

0.33 

319 

0.64 

0.11 

0.61 

0.05 

45 

H-3L 

0.36 

5.48 

0.40 

378 

1.07 

0.10 

0.73 

0.11 

44 

H-4L 

0.36 

7.37 

0.53 

509 

1.15 

0.15 

1.09 

0.08 

41 

H-5L 

0.36 

6.94 

0.83 

479 

0.72 

0.10 

1.22 

0.82 

32 

H-IM 

0.73 

2.96 

0.21 

204 

0.74 

0.06 

0.62 

0.18 

66 

H-2M 

0.73 

5.04 

0.28 

348 

1.23 

0.37 

1.16 

0.80 

63 

H-3M 

0.73 

2.07 

0.19 

143 

0.27 

0.03 

0.24 

0.05 

49 

H-4M 

0.73 

3.22 

0.19 

222 

0.75 

0.17 

0.69 

0.31 

47 

H-5M 

0.73 

4.30 

0.28 

297 

0.58 

0.10 

0.34 

0.02 

43 

H-IH 

1.38 

3.01 

0.19 

208 

2.34** 

0.24 

j  47** 

0.54 

85 

H-2H 

1.38 

3.64 

0.26 

251 

0.86 

0.07 

0.53 

0.09 

80 

H-IZ 

0 

8.72 

0.96 

602 

2.41* 

2.03 

0.86 

0.15 

62 

H-2Z 

0 

9.93 

0.45 

685 

1.80 

0.43 

1.98 

0.31 

59 

H-3Z 

0 

11.64 

0.73 

803 

3.10 

1.04 

2.95 

0.59 

30 

Table  5-2:  Vertical  Diffusivity  Results  for  Field  Study 


Run# 

P 

del  u 

(cm/s) 

Re-d 

D-MOM 

(cm''2/s) 

delD 

(cm''2/s) 

D-EF 

(cm''2/s) 

delD 

(cm^2/s) 

Depth 

(cm) 

V-IL 

0.36 

5.47 

0.58 

377 

0.58 

0.20 

0.51 

0.06 

45 

V-2L 

0.36 

6.84 

0.56 

472 

1.28 

0.19 

0.92 

0.16 

43 

V-3L 

0.36 

6.35 

0.47 

438 

1.37 

0.14 

1.37 

0.10 

44 

V-4L 

0.36 

7.57 

0.44 

522 

1.12 

0.08 

1.25 

0.50 

38 

V-IM 

0.73 

3.25 

0.20 

224 

0.78 

0.05 

0.64 

0.10 

66 

V-2M 

0.73 

3.46 

0.15 

239 

0.53 

0.04 

0.39 

0.09 

61 

V-3M 

0.73 

2.33 

0.23 

161 

0.36 

0.04 

0.18 

0.04 

49 

V-4M 

0.73 

3.76 

0.14 

259 

0.58 

0.03 

0.53 

0.04 

46 

V-IH 

1.38 

3.96 

0.33 

273 

0.71 

0.06 

0.52 

0.05 

77 

V-IZ 

0 

9.52 

0.89 

657 

0.69 

0.07 

0.70 

0.10 

33 

V-2Z 

0 

10.25 

0.88 

707 

1.19 

0.49 

1.07 

0.55 

31 

V-3Z 

0 

13.10 

0.65 

904 

1.84 

0.15 

1.96 

0.24 

29 

Notes: 

*  -  Diffusivity  calculated  using  the  Method  of  Moments  for  profile  H-IZ  is  anomalously  high  due  to  the 
coordinates  of  the  profile  points  and  the  extrapolation  technique  used.  This  diffusivity  value  was  not  used 
in  trend  analysis,  although  it  is  plotted  in  figures. 

**  -  Diffusivities  calculated  for  profile  H-IH  are  also  thought  to  be  anomalous  based  on  comparison  to 
other  giniilar  flow-vegetation  combinations.  These  values  were  also  not  used  in  trend  analysis,  but  again  are 
plotted  in  applicable  figures. 
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5.2  Vegetated  versus  Non-Vegetated  Diffusivity 


One  of  the  primary  focuses  of  this  research  was  examination  of  the  effects  of  vegetation 
on  diffusivity.  Comparison  of  diffiisivities  observed  for  vegetated  and  non-vegetated 
profiles  provides  a  direct  means  of  examining  this  effect.  This  section  describes  the 
methods  used  to  make  this  comparison  and  the  conclusions  reached. 

5.2.1  Non-Vegetated  Diffusivity 

Turbulent  diffiisivities  are  known  to  depend  on  turbulence  intensity  and  some 
characteristic  turbulent  length  scale  (Fischer,  et  al,  1979).  Flow  in  non-vegetated  open 
channels  has  been  widely  studied,  and  several  empirical  relationships  relating  diffusivity 
to  channel  flow  and  geometry  have  been  developed.  Fischer,  et  al  (1979,  p.  1 12)  give 
Equation  5-1  as  a  means  of  estimating  transverse  diffiisivities  in  meandering,  irregular¬ 
sided  channels.  This  empirical  relation  was  developed  using  data  fi-om  seven  separate 
studies  and  has  a  stated  accuracy  of  ±50%.  The  constant  of  proportionality  in  Equation 
5-1  is  reduced  to  0.15  for  smooth,  rectangular  channels  and  by  an  order  of  magnitude  (to 
0.06)  for  vertical  diffiisivities  in  natural  channels  (Fischer,  et  al,  1979). 

s,  =  0.6hu  *  Equation  5-1 

(h  is  the  mean  depth  of  the  chaimel  and  u*  is  tiie  shear  velocity  at  the  bed.) 

Elder  (1959)  derived  a  similar  expression  for  vertical  diffiisivities  in  turbulent  shear  flow 
over  a  wide,  inclined  plane.  Equation  5-2  gives  this  expression.  This  relation  has  been 
widely  proven  and  accepted,  and  matches  well  with  the  expected  order  of  magnitude 
reduction  from  the  Equation  5-1  expectation  as  discussed  above. 

=  0.067 hu  *  Equation  5-2 

Comparison  of  measured  diffiisivities  with  these  empirically  expected  results  was  the  first 
analysis  undertaken.  Shear  velocity  was  not  measured  during  the  experiment  and  so  had 


to  be  estimated.  For  the  seven  studies  used  to  develop  Equation  5-1  the  mean  shear 
velocity  was  7%  of  mean  flow  (Fischer,  et  al,  1979),  so  this  was  used  as  a  reasonable 
estimate  of  shear  flow  for  the  non-vegetated  profiles.  Using  this  estimated  shear  flow  and 
Equations  5-1  and  5-2,  the  slopes  of  diffiisivity  plotted  versus  u-h  were  expected  to  be 
0.04  ±  0.02  for  transverse  difihisivity  and  0.005  for  vertical  diffiisivity. 

Plotting  D  vs.  u-h  for  unvegetated  profiles  and  then  applying  a  linear  least  squares  fit 
(with  forcing  through  (0,0))  resulted  in  slopes  of  0.004  for  both  transverse  and  vertical 
diffiisivity.  This  is  within  20%  of  the  expected  value  for  vertical  diffiisivity,  but  is  an 
order  of  magnitude  less  than  that  expected  for  transverse  diffiisivity.  This  discrepancy  is 
believed  to  be  due  to  the  scales  involved  in  the  current  study.  Because  the  short  temporal 
and  spatial  scales  used  were  too  small  to  resolve  the  longer-term  variations  in  currents 
and  turbulence  structure  which  would  tend  to  increase  diffiisivity,  measured  transverse 
dififiisivities  may  have  been  too  small.  Transverse  and  vertical  scales  were  essentially  the 
same  in  this  study,  since  mean  depth  was  on  the  order  of  0.5  meters  and  the  transverse 
plume  dimensions  were  roughly  the  same  size.  This  is  a  marked  difference  firom  the 
roughly  100:1  widfii  to  depth  ratio  for  the  studies  used  to  develop  Equation  5-1  (Fischer, 
et  al,  1979),  and  this  may  accoimt  for  die  lower  than  expected  transverse  diffiisivities. 
This  argument  is  further  strengthened  by  reference  to  another  relation  firom  Fischer  (1979, 
p.77)  showing  how  observed  diffiisivities  change  with  scale.  Diffiisivities  (in  cm  /s)  are 
shown  to  be  related  to  dye  cloud  size  (in  cm)  raised  to  the  4/3  power  by  a  constant  of 
proportionality  between  0.002  and  0.01  (Fischer,  et  al,  1979).  Applying  this  relation  to  a 
0.5  meter  dye  cloud  provides  estimated  diffiisivities  from  0.4  to  1 .9  cm^/s.  Again  this 
matches  well  with  measured  vertical  diffiisivities,  aldiough  measured  transverse 
diffiisivities  (ranging  between  1  and  3  cm^/s)  are  now  slightly  higher  than  diese  new 
expected  values.  These  differences  are  minor  and  can  easily  be  accounted  for  by  small 
variations  in  dye  cloud  size,  however.  Measured  transverse  diffiisivities  therefore  seem 
to  be  reasonable  based  on  expected  results  for  the  small  scales  employed. 
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5.2.2  Vegetated  Diffusivity 


The  technique  of  relating  D  to  u-h  was  also  applied  to  vegetated  profiles  to  examine  the 
influence  of  vegetation  on  diffusivity.  With  the  addition  of  vegetation,  both  the  length 
scale  and  turbulence  intensity  which  control  diffusion  are  expected  to  change  firom  the 
non-vegetated  case.  The  length  scale  can  be  expected  to  change  firom  the  chaimel  depth 
to  somewhere  on  the  order  of  mean  stem  diameter  or  mean  stem  spacing  for  the 
vegetation  (Ackerman  and  Okubo,  1993,  Anderson  and  Charters,  1982,  Nepf,  et  al,  1997). 
If  nothing  else  changes,  diffusivities  in  vegetation  would  then  be  expected  to  be  reduced 
relative  to  non-vegetated  diffusivities  based  on  this  change  in  length  scale.  There  are 
other  changes,  however,  most  notably  in  the  turbulence  intensity.  In  the  non-vegetated 
case,  bed-generated  turbulence  dominates  and  shear  velocity  at  the  bed  governs 
turbulence  intensity.  In  the  vegetated  case  this  bed-generated  turbulence  is  supplemented 
by  stem-generated  turbulence,  however,  so  it  is  no  longer  valid  to  use  shear  velocity  as  a 
measure  of  turbulence  intensity  or  as  a  scaling  factor.  Both  bed-  and  stem-generated 
turbulence  can  be  expected  to  scale  on  the  mean  flow,  however,  so  mean  flow  can  be 
used  instead.  Turbulence  intensity  is  then  given  by  some  constant  times  the  mean  flow, 
and  diffusivity  is  given  by  another  constant  of  proportionality  multiplied  by  this 
turbulence  intensity  and  the  turbulent  length  scale.  The  end  result  is  Equations  5-3  and 
5-4,  where  the  two  constants  of  proportionality  have  been  combined  into  a  single  A  and  A ' 
depending  on  whether  mean  depth  (h)  or  mean  stem  diameter  {d)  is  used  as  the 
representative  turbulent  length  scale.  As  noted  above,  previous  studies  have  shown  d  to 
be  the  correct  length  scale  within  vegetation,  so  Equation  5-4  is  expected  to  be  more 
correct. 


D=A'U'h  Equation  5-3 

D—A  ''U'd  Equation  5-4 

For  more  direct  comparison  to  the  non-vegetated  case.  Equation  5-3  was  first  used  to 
analyze  data  from  the  vegetated  profiles.  D  was  plotted  against  u  h  for  all  vegetation 


densities  separately  and  in  combination,  with  a  1 
forcing  through  (0,0))  used  to  iSnd  the  best  fit  va 
Figure  5-1  and  Figure  5-2,  with  results  tabulated 
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►  Figure  5-1 :  D  vs.  u-h  for  Transverse  Profiles 
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Figure  5-2:  D  vs.  u*h  for  Vertical  Profiles 


ir  least  squares  technique  (wifii 
for  A  in  each  case.  Data  is  plotted 


Table  5-3:  Best  Fit 

Slopes  for  D  vs.  u 

•h  Plots 

Profiles  Used 

Transverse  Slope 

Vertical  Slope 

Unvegetated 

0.0040 

0.0038 

Low  Density 

0.0035 

0.0038 

Medium  Density 

0.0036 

0.0028 

High  Density 

0.0046 

0.0020 

All  Vegetated 

0.0035 

0.0033 

This  analysis  shows  that  vegetated  and  non-vegetated  diffiisivities  are  very  comparable. 

In  general,  the  computed  slopes  are  smaller  for  vegetated  runs  than  for  non-vegetated 
runs,  meaning  that  diffusivity  in  the  vegetated  areas  would  be  less  than  for  an  open 
channel.  This  decrease  for  vegetated  conditions  is  quite  small,  however,  and  one  slope 
actually  increased  while  another  remained  unchanged.  Therefore,  based  on  the  small 
differences  measured  and  the  imcertainty  in  the  data  gathered,  it  is  difiBcult  to  say  whether 
any  significant  difference  exists  between  vegetated  and  unvegetated  diffiisivities  based  on 
this  analysis.  This  result  is  based  on  using  water  depdi  as  a  scaling  factor.  Similar 
analysis  techniques  using  different  characteristic  length  scales  were  used  next  to  further 
examine  any  change  in  diffiisivities  due  to  vegetation. 

The  preceding  analysis  scaled  vegetated  diffusivity  with  mean  water  depth  for  ease  of 
comparison  with  the  unvegetated  cases,  but  as  noted  earlier  vegetated  diffusivity  would 
be  more  aptly  scaled  by  mean  stem  diameter  or  mean  spacing  between  the  stems.  Nepf, 
et  al  (1997),  using  vertical  wooden  dowels  to  simulate  vegetation  stems,  found  that  the 
characteristic  length  scales  for  turbulence  within  the  dowels  (stems)  was  on  the  order  of 
the  stem  diameter.  Mean  stem  diameter  will  be  examined  first  as  a  scaling  factor  in  this 
study,  with  analysis  using  mean  stem  spacing  to  follow.  Since  d  in  Equation  5-4  is  now 
a  constant,  analysis  based  on  stem  diameter  reduces  to  evaluating  diffusivity  versus  mean 
flow.  Mean  flow  was  converted  to  a  Reynolds  number  based  on  mean  stem  diameter  and 
least  squares  analyses  were  conducted  as  before.  The  resulting  data  is  plotted  in  Figure  5- 
3  and  Figure  5-4  and  tabulated  in  Table  5-4.  Uncertainty  values  given  represent  the  95% 
confidence  level,  and  non-vegetated  data  is  plotted  for  comparative  purposes. 


Diffusivity  (cm^2/s)  ^  Diffusivity  (cm'"2/s) 


i  5-3:  Transverse  and  Vertical  Diffusivities  vs.  Reynolds  Number 


o  D-  Low 
■  D  -  Med 
A  D- 

•  D-  Zero 


Reynolds  Number  (based  on  stem  diameter) 


Figure  5-4:  Diffusivities  vs.  Reynolds  Number  by  Vegetation  Density 
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Table  5-4:  Best  Fit  Slopes  for  D  vs.  Re  Plots 


Transverse 

Vertical  | 

Data  Points  Used 

Slope  of  D  vs.  Re 

Uncertainty 

Slope  of  D  vs.  Re 

Uncertainty 

All  vegetated  data 

0.0023 

0.0005 

0.0023 

0.0003 

Low  density  points 

0.0020 

0.0003 

0.0024 

0.0004 

Medixun  density  points 

0.0028 

0.0006 

0.0023 

0.0005 

High  density  points 

0.0028 

0.0013 

0.0023 

0.0007 

As  can  be  seen,  best  fit  slopes  for  both  transverse  and  vertical  dififusivity  were  essentially 
equal  to  each  other  over  all  vegetation  densities.  Transverse  diffusivities  at  medium  and 
high  vegetation  densities  seem  to  be  slightly  higher  than  other  values,  but  their 
uncertainty  is  also  higher.  From  this  data,  it  appears  that  the  vegetation  density  has  little 
effect  on  measured  diffiisivities. 


The  mean  stem  spacing  was  also  examined  as  a  turbulence  length  scale.  Stem  spacing 
should  be  inversely  proportional  to  vegetation  density,  so  diffusivity  was  plotted  versus 
u/P  to  examine  dependence.  Figure  5-5  provides  the  results.  Slopes  of  the  best  fit  D  vs. 
u/P  lines  shown  are  0.05, 0.13,  and  0.24  for  low,  medium,  and  high  stem  densities.  These 
differences  may  be  significant,  but  the  data  are  perhaps  more  indicative  of  the  relation 
between  stem  density  and  current  velocity.  Points  for  each  stem  density  are  clearly 
clustered  in  specific  u/P  ranges,  showing  that  stem  density  controls  flow. 


♦  Low 
□  Medium 
A  High 

"  "  ■  Linear  (Low) 

Linear  (Medium) 
“  — Linear  (High) 


Figure  5-5:  Diffusivity  vs.  u/P  for  Vegetated  Profiles 
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5.3  Transverse  versus  Vertical  Diffusivity 


Figure  5-3,  Figure  5-4,  and  Table  5-4  also  show  the  close  agreement  between  transverse 
and  vertical  diffusivities.  They  both  appear  to  be  directly  related  to  ambient  flow,  and 
their  values  at  similar  flows  seem  to  be  comparable  across  all  vegetation  densities.  As 
can  be  seen,  there  is  very  little  difference  in  measured  transverse  and  vertical  diffusivities, 
particularly  at  low  Reynolds  numbers.  Up  to  a  Re  of  roughly  500  the  values  fall  ri^it  on 
top  of  one  another,  although  at  hi^er  flows  some  divergence  between  the  two  values 
may  be  evident.  Data  in  this  region  is  much  more  sparse  and  uncertain,  however,  so  this 
trend  may  or  may  not  be  significant.  The  slopes  computed  using  data  for  all  vegetated 
profiles  were  the  same  for  transverse  and  vertical  diffusivity,  and  slopes  computed  for 
individual  vegetation  densities  all  agreed  within  the  95%  confidence  levels  given. 

This  close  agreement  between  the  transverse  and  vertical  diffusivities  was  not  necessarily 
expected,  but  is  also  not  necessarily  too  surprising  either.  A  compilation  of  turbulent 
diffusion  studies  made  in  natural  channels  (Fischer,  et  al,  1979)  shows  that,  in  general, 
vertical  diffusivities  are  roughly  an  order  of  magnitude  less  than  transverse  diffusivities. 
This  difference  is  due  to  the  different  length  scales  of  turbulence  in  the  transverse  and 
vertical  directions  -  in  natural  channels  (with  no  plants)  there  is  a  wide  difference  in  these 
two  length  scales.  Laboratory  studies  (Nepf,  et  al,  1997)  using  vertical  wooden  dowels  to 
simulate  plants  have  shown  a  marked  reduction  in  this  length  scale  difference,  and  hence 
in  diffusivity.  The  dowels  (plant  stems)  act  to  break  up  the  turbulent  eddies  oriented  in 
the  transverse  direction,  with  the  characteristic  length  scale  being  reduced  to  on  the  order 
of  the  spacing  between  the  dowels  (Nepf,  et  al,  1997).  Transverse  and  vertical 
dififiisivities  differ  by  roughly  a  factor  of  four  in  these  laboratory  studies,  a  sigmficant 
reduction  from  the  factor  of  ten  found  in  natural,  plant-firee,  systems  (Fischer,  et  al,  1979). 

In  this  study,  the  plant  stems  should  act  much  like  the  dowels  to  reduce  tihie  transverse 
turbulent  length  scale.  Unlike  the  dowels,  however,  in  the  natural  setting  the  stems  (and 
leaves)  grow  in  many  directions  centered  aroxmd  the  vertical.  This  spreading  of  the  stem 
matrix  should  produce  wakes  with  a  wide  variety  of  orientations,  therefore  reducing  the 
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2-D  bias  of  the  flume  experiment  and  leading  to  reduced  anisotropy  in  the  transverse  and 
vertical  directions.  Additionally,  the  natural  stems  are  quite  flexible  and  move  under  the 
influence  of  the  wind  and  the  current.  This  should  enhance  the  stirring  effect  and  lead  to 
more  similar  transverse  and  vertical  turbulence  scales.  These  reasons,  coupled  with  the 
scale  discussion  presented  earlier,  makes  the  close  agreement  between  the  transverse  and 
vertical  diffusivity  values  seems  reasonable. 

5.4  Vegetation  Density  Effects 

As  seen  from  Figure  5-4  and  Table  4-3,  there  is  very  little  difference  in  the  dififusivities 
observed  at  different  vegetation  densities.  One  must  be  careful  in  making  the  statement 
that  difiEusivity  is  ^  dependent  on  vegetation  density  from  this  data,  however,  since  there 
are  very  few  cases  where  different  vegetation  densities  are  sampled  at  similar  or 
overlapping  Reynolds  numbers. 

These  results  differ  somewhat  from  those  obtained  m  the  laboratory  study  discussed 
earlier  (Nepf,  et  al,  1997).  Figure  5-6  presents  data  from  this  laboratory  study,  and  shows 
quite  clearly  that  increasing  the  dowel  (stem)  density  increases  diffusivity  for  a  given 
Reynolds  number  in  that  experiment.  Diffusivity  appears  to  increase  linearly  with 
Reynolds  number  from  these  results,  similar  to  the  conclusions  reached  in  the  current 
study.  The  increase  in  diffusivity  for  a  given  increase  in  dowel  density  is  also  seen  to  fall 
as  the  dowel  density  increases  due  to  shading  effects.  Again,  this  increase  in  diffusivity 
with  increasing  stem  density  was  not  seen  in  the  field  study.  (For  reference.  High 
vegetation  density  in  this  field  study  was  1.4%,  Medium  was  0.7%,  and  Low  was  0.35%). 


Figure  5-6;  Transverse  Diffusivity  versus  Reynolds  Number  -  Laboratory  Study 
(data  from  Nepf,  et  al,  1997) 


There  are  several  explanations  as  to  why  results  from  this  study  do  not  show  the 
dependence  of  diffusivity  on  vegetation  density  evident  in  the  laboratory  study.  The 
laboratory  study  was  able  to  measure  diffusivity  for  different  stem  densities  at  similar 
Reynolds  numbers,  but  in  the  field  this  was  not  done  to  any  great  extent.  The  added  flow 
restriction  imposed  by  the  additional  vegetation  was  sufficient  to  re-route  flow  around  the 
vegetation  patch,  so  any  diffusivily-enhancing  effect  of  the  increased  stem  density  tended 
to  be  counterbalanced  by  the  concurrent  flow  reduction  within  the  patch.  These  two 
effects  essentially  cancel  each  other  out,  and  little  diffusivity  dependence  on  stem  density 
is  noted.  Vegetation  density  seems  to  control  flow,  and  this  flow  then  determines 
diffusivity. 

Alternately,  it  could  be  possible  tiiat  the  flexibility  and  the  more  random  orientation  of  the 
natural  stems  (relative  to  the  dowels)  magnifies  the  shading  effect  seen  in  the  laboratory 
data  at  increasing  dowel  density.  Since  only  a  few  stems  (and  leaves)  would  then  be 
necessary  to  effectively  cover  the  full  cross  section  for  flow,  it  may  be  that  there  indeed  is 
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not  much  diffiisivity  enhancement  as  vegetation  density  is  increased.  If  this  were  true,  the 
presence  of  only  a  few  plants  would  increase  the  dififusivity  significantly,  but  addition  of 
more  plants  would  have  little  effect. 

Another  possible  explanation  for  the  seeming  independence  on  vegetation  density  is  that, 
at  least  for  this  site,  other  factors  dominate.  The  site  was  chosen  so  as  to  minimize 
turbulence  induced  by  means  other  than  the  vegetation,  but  assessment  of  these  other 
turbulence  sources  may  have  underestimated  their  contribution.  Although  the  region 
around  the  vegetation  patch  did  have  a  fairly  smooth  bed,  some  roughness  and 
xmevenness  was  still  present,  so  bed-generated  turbulence  may  have  been  more  important 
than  originally  thought.  In  addition,  the  patch  lay  rou^y  25-30  meters  downstream  of  a 
major  channel  bend  on  both  incoming  and  outgoing  tides.  Tinbulent  mixing  coefficients 
have  been  measured  at  2-25  times  their  normal  value  downstream  of  bends  in  natural 
channels  (Fischer,  et  al,  1979),  so  perhaps  significant  residual  turbulence  from  these 
bends  still  existed  at  the  sampling  point.  Contributing  to  the  effects  of  both  these 
concerns,  only  roughly  1.5  meters  of  vegetation  existed  as  a  zone  of  transition  between 
die  free  stream  and  the  dye  injection  point.  This  distance  was  the  maximum  available  at 
this  particular  vegetation  patch,  but  was  perhaps  not  sufficient  to  fully  smooth  the  flow 
before  beginning  the  experiment.  These  three  effects;  turbulence  contributed  by  the  bed, 
the  natural  turbulence  enhancement  due  to  the  bends  in  the  channel,  and  insufficient 
smoothing  of  the  flow  within  the  patch;  when  taken  together  may  have  overwhelmed  the 
contribution  of  the  differing  plant  densities,  making  it  seem  that  diffiisivity  is  independent 
of  vegetation  density. 

To  examine  the  dependence  of  diffiisivity  on  vegetation  density  a  little  more  closely,  the 
non-dimensional  parameter  D/u-d  was  plotted  versus  percent  stem  coverage.  Figures  5-7 
and  5-8  show  these  relationships,  with  the  open  symbols  being  the  individual  data  points 
and  the  solid  symbols  being  the  mean  value  at  each  vegetation  density.  Although  they 
may  not  be  significant  due  to  the  uncertainties  involved,  the  results  from  this  examination 
seem  to  support  the  conclusions  reached  earlier.  D/u-d  is  essentially  constant  across  all 
vegetation  densities  for  transverse  diffiisivity,  although  there  may  be  a  slight  reduction  in 


the  vegetated  versus  the  unvegetated  cases.  For  vertical  diffusivity,  there  is  somewhat  of 
an  increase  in  D/u-d  when  going  from  unvegetated  to  low  density,  but  then  the  value 
remains  constant  as  vegetation  density  is  increased.  These  results  support  the 
enhancement  of  vertical  difiEusivity  by  addition  of  a  few  stems,  followed  by  little  change 
as  more  stems  are  added,  and  the  slight  reduction  in  horizontal  difEiisivity  found  in  the 
previous  analysis.  Values  of  D/u-d  between  transverse  and  vertical  dififusivities  also  were 
quite  comparable  to  one  another,  again  confirming  earlier  results. 


Figure  5-7:  D/ud  vs.  Percent  Stem  Coverage  -  Transverse  Diffusivity 


Figure  5-8:  D/ud  vs.  Percent  Stem  Coverage  -  Vertical  Diffusivity 
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As  a  further  examination  of  the  influence  of  vegetation  density  on  diffusivity,  D  was  also 
plotted  as  a  function  of  P  Re  for  vegetated  profiles.  These  relationships  for  horizontal 
and  vertical  diffusivity  are  shown  in  Figures  5-9  and  5-10,  respectively,  with  the  solid  line 
being  a  least  squares  fit  to  the  data.  Although  there  is  again  much  uncertainty,  it  appears 
that  diffusivity  values  for  PRe  greater  than  about  200  begin  to  fall  off  from  the  best  fit 
line.  It  appears  the  initial  slope  is  roughly  0.005  for  both  cases,  but  that  slopes  obtained 
using  only  the  higher  P  Re  points  would  be  much  lower.  This  observation  is  likely  due  to 
the  shading  effects  discussed  earlier,  and  could  come  about  as  either  vegetation  density  or 
flow  is  increased  sufficiently.  Increasing  the  stem  density  obviously  increases  the 
chances  of  shading  as  there  are  more  stems  to  interact  with.  Increasing  the  flow  causes  an 
increase  in  turbulence  and  vortex  shedding,  and  the  enlarged  vortices  have  a  higher 
chance  of  interacting  with  other  stems.  P  Re  =  200  appears  to  be  the  breaking  point  for 
reduction  of  diffusivity,  although  further  study  would  be  necessary  to  say  this  with  any 
certainty. 


Diffusivity  (cm'^2/s) 


Figure  5-9:  D  vs.  P*Re  for  Transverse  Profiles 
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5.5  Other  Comparisons 


Difiusivity  values  obtained  for  vegetated  profiles  do  compare  well  with  other  field  and 
laboratory  experiments  done  in  vegetation.  A  mixing  study  done  at  low  flow  in  eelgrass 
beds  found  a  mixing  coefficient  of  roughly  0.5  +  0.4  cm^/sec  at  cmnrent  velocities  within 
the  bed  of  2  cm/sec  (Worcester,  1995).  This  flow  is  equivalent  to  a  Reynolds  number  of 
170  based  on  mean  stem  diameter,  and  the  diffusivity  value  above  matches  well  with 
those  presented  in  Figure  5-3.  The  dye  patch  scales  used  in  that  study  were  0. 1  m  to 
several  meters  and  little  difference  was  noted  between  vegetated  and  non-vegetated 
diffusivities  (Worcester,  1995),  again  matching  well  with  this  experiment.  Measured 
values  from  the  current  field  experiment  also  match  well  with  the  laboratory  values 
obtained  by  Nepf,  et  al,  in  their  dowel-flume  experiment  (see  Figure  5-6).  Although  the 
uncertainty  in  the  data  from  this  field  experiment  is  much  greater,  the  magnitudes  of  the 
individual  points  agree  fairly  well,  at  least  for  the  vegetated  profiles.  The  diffusivity 
found  in  this  experiment  for  the  unvegetated  profiles  is  much  greater  than  that  found  in 
the  lab,  but  there  are  important  differences  in  the  two  studies  as  well.  The  flume  was  a 
limited  size  in  width  and  hei^t,  therefore  the  eddy-scale  is  much  smaller  than  in  the 
field.  Also,  the  bed  in  the  field  is  much  less  smooth  than  in  the  laboratory,  so  additional 
contributions  from  bed-generated  turbulence  can  be  expected. 

5.6  Summary  and  Recommendations 

In  summary,  the  results  obtained  for  diffusivity  appear  to  be  reasonable  based  on  intuition 
and  previously  released  results.  The  goal  of  this  study  was  to  examine  the  effects  of 
vegetation  density  on  diffusion,  but  the  results  showed  little  such  dependence. 
Diffusivities  were  found  to  be  much  more  dependent  on  flow  than  on  vegetation  density. 
Increased  vegetation  density  does  cause  a  reduction  in  flow,  thereby  reducing  diffusivity 
as  a  secondary  effect.  Diffusivity  withm  vegetated  regions  is  therefore  expected  to  be 
reduced  relative  to  non-vegetated  regions.  Transverse  and  vertical  diffusivities  were  also 
found  to  be  isotropic  across  all  ranges  of  vegetation  density,  a  departure  from  earlier 
published  results.  This  isotropy  is  accomplished  even  at  low  stem  densities  and  is 
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believed  to  be  due  to  the  nearly  random  orientation  of  the  stems.  Laboratory  studies  and 
models,  such  as  those  by  Nepf,  et  al,  (1997),  need  to  be  revised  to  account  for  this  result. 
Additionally,  techniques  used  in  this  study  were  partially  ineffective  in  that  sampling 
times  were  too  small  to  resolve  longer  time  scale  variations  in  diffusion.  To  address  these 
issues,  suggestions  for  improvements  in  further  studies  include: 

•  Future  studies  should  be  done  at  a  site  which  provides  a  longer  time  period  for 
sampling.  This  could  be  done  by  finding  a  tidal  marsh  where  the  water  level  does 
not  rise  and  fall  as  quickly  as  it  did  at  this  location,  or  perhaps  even  by  locating  an 
alternate  spot  in  Great  Sippewissett  Marsh  closer  to  the  inlet  channel.  These 
options  are  likely  to  offer  only  limited  improvement,  however,  so  a  better  option 
may  be  to  conduct  a  study  in  a  vegetated  creek  or  river  channel.  This  would 
provide  consistent  study  conditions  over  the  period  necessary  to  resolve  longer- 
term  variability. 

•  Alternately,  a  different  sampling  technique  could  be  used  to  effectively  gather 
more  data  at  each  point  over  the  same  available  time  period.  This  would  entail 
sampling  from  several  locations  concurrently,  significantly  increasing  the  time, 
effort,  and  equipment  required  for  sampling.  One  possible  approach  is  to 
simultaneously  collect  water  samples  fi'om  a  number  of  points  using  a  number  of 
small  pumps,  and  then  analyzing  these  water  samples  for  fluorescence  either  in 
the  lab  or  in  the  field  in  discrete  sample  mode.  This  technique  would  provide  a 
longer-term  average  fluorescence  value,  but  all  time  dependence  would  be  lost. 
Another  possibility  is  to  use  several  fluorometers  simultaneously,  but  again  this 
significantly  expands  the  effort  and  expenditure  required.  Fluorescence  readings 
in  the  field  would  also  be  much  more  preferable,  since  transporting  large  water 
samples  would  be  quite  difBcult. 

•  Additionally,  some  provision  needs  to  be  made  for  ensuring  comparable  flows 
through  the  study  site  at  different  vegetation  densities.  In  this  study,  even 
relatively  sparse  vegetation  coverage  was  sufficient  to  reroute  flow  arovmd  the 
vegetation  patch,  so  that  increased  diffusivity  due  to  increased  stem  density  was 


masked  by  the  concurrent  flow  reduction.  This  made  it  difficult  to  assess  the 
effect  of  the  vegetation  coverage,  since  very  few  points  at  similar  flows  but 
differing  vegetation  densities  were  gathered.  To  adequately  measure  the  true 
effects  of  vegetation  density  on  diffiisivity,  a  site  where  the  water  is  constrained  to 
flow  throu^  the  vegetated  area  would  be  much  better.  Again,  a  vegetated  creek 
or  river  channel  would  work  well,  particularly  one  vegetated  from  bank  to  bank. 

•  Finally,  more  and  better  instrumentation  and  recording  is  probably  warranted  in 
future  studies.  Currents  obtained  in  this  study  were  adequate  for  the  calculations 
which  were  done,  but  were  quite  crudely  obtained.  Use  of  a  continuously 
operating,  low-flow  current  measurement  device,  such  as  a  BASS  (Benthic 
Acoustic  Stress  Sensor)  instrument,  would  allow  better  resolution  of  current 
variability  and  magnitude  over  the  coiurse  of  a  profile.  This  would  enable  better 
and  more  detailed  analysis  of  the  results,  and  may  help  in  correlation  between 
current  variability  and  fluorescence  fluctuations.  Provisions  should  also  be  made 
for  more  carefully  recording  water  depth  during  the  course  of  the  experiment, 
perhaps  by  the  use  of  a  floating  tidal  height  gauge.  Finally,  more  careful 
measurement  of  the  vertical  coordinates  involved  (such  as  dye  injection  hei^t 
and  sample  port  location)  would  allow  more  full  use  of  the  data  obtained  for 
proper  curve  fitting.  Since  vertical  coordinates  were  not  well  resolved  in  the 
current  study,  use  of  two-dimensional  contouring  and  analysis  of  the  leading 
constant  in  the  one-dimensional  exponential  curve  fitting  were  not  done. 

Addition  of  these  capabilities  would  greatly  enhance  future  studies. 
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Appendix  A  -  Detailed  Results 


This  appendix  presents  more  detailed  data  and  results  for  the  study.  Detailed  tables 
showing  data  points  used  and  intermediate  results  for  each  profile  are  included  first, 
followed  by  plots  of  the  individual  fluorescence  profiles. 

Each  profile  presented  here  is  identified  by  a  three  character  ‘Run  Number’  designation. 
This  designation  begins  with  ‘H’  for  horizontal  profiles  and  with  ‘V’  for  vertical  profiles. 
The  second  character  is  an  incremental  counter,  starting  at  ‘1’  for  horizontal  and  vertical 
profiles  at  each  plant  density  and  increasing  for  each  successive  horizontal  or  vertical 
profile.  The  third  character  identifies  the  vegetation  density  associated  with  the  profile, 
and  is  either  an  ‘L’,  ’M’,  ‘H’,  or  ‘Z’  depending  on  whether  the  profile  is  for  low, 
medium,  high,  or  zero  plant  density.  All  profiles  were  analyzed  using  the  techniques 
described  in  tiie  thesis,  with  the  exception  of  the  following  variations: 

•  Profile  H-IZ  -  The  mean  fluorescence  values  for  the  first  and  second  data 
point  in  this  profile  were  very  nearly  equal.  These  fluorescence  readings  also 
had  uncertainties  quite  large  relative  to  their  actual  value.  These  two  things 
acting  together  made  using  these  points  to  extrapolate  to  zero  fluorescence 
impossible  -  fluorescence  never  reached  zero  when  using  certain  combinations 
of  these  points.  The  combinations  which  caused  this  problem  were  not 
included  in  the  analysis  of  this  profile. 

•  Profiles  V-IZ,  V-2Z,  and  V-3Z  -  These  profiles  were  all  taken  in  relatively 
shallow  water  and  with  the  dye  source  at  or  above  the  water  surface.  Because 
the  dye  was  distributed  only  in  the  water,  however,  only  half  of  the  plume 
described  by  the  equations  used  actually  developed  for  these  profiles.  In  these 
cases,  a  sample  was  taken  as  close  to  the  water’s  surface  as  possible  (roughly 
2  cm)  without  drawing  air  into  the  sample  pump.  The  water  siuface  was  then 
set  as  a  mirror  at  zero  elevation  and  the  measured  points  were  reflected  in  the 


mirror  to  give  a  full  profile.  The  measured  point  closest  to  the  water  was 
given  a  coordinate  of  2  cm  displacement  firom  the  mirror,  with  other 
coordinates  adjusted  similarly.  The  resulting  full  profiles  were  then  analyzed 
with  the  standard  methods  described  above.  Transformed  coordinates  were 
converted  back  to  their  original  values  for  plotting,  with  only  the  lower  half  of 
the  profile  (the  portion  actually  in  the  water)  being  plotted. 


Table  A-1:  Profile  Data  for  Low  Plant  Density 
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Table  A-3:  Profile  Data  for  High  Plant  Density 
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Table  A-5:  Diffusivity  Results  for  Low  Plant  Density 
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Table  A-6:  Diffusivity  Results  for  Medium  Plant  Density 
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Table  A-8:  Diffusivity  Results  for  Zero  Plant  Density 
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Figure  A-2:  Vertical  Diffusivity  Results  -  Low  Plant  Density 
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Figure  A-3:  Horizontal  Diffusivity  Results  -  Medium  Plant  Density 
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80 


Diffusivity  ( 


Figure  A-7:  Horizontal  Diffusivity  Results  -  Zero  Plant  Density 


Figure  A-8:  Vertical  Diffusivity  Results  -  Zero  Plant  Density 
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Fluorescence  (ppb)  ^  |  Fluorescence  (ppb) 


A-9:  Profile  V-IL  -  August  6th,  Set  B,  Points  7-9 
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Figure  A-10:  Profile  V-2L  -  August  6th,  Set  C,  Points  1-4 
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Fluporescence  (ppb)  n  \  Fluo 


A-ll:  Profile  H-IL  -  August  6th,  Set  C,  Points  6-11 


Figure  A-12:  Profile  H-2L  -  August  7th,  Set  A,  Points  1-4 
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Fluorescence  (ppb)  ^  Fluorescence  (ppb) 


A-13:  Profile  H-3L  -  August  7th,  Set  A,  Points  5-10 
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Figure  A-14:  Profile  V-3L  -  August  7th,  Set  A,  Points  11-14 
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Fluorescence  (ppb)  ^  |  Fluo 


!  A-15:  Profile  H-4L  -  August  7th,  Set  B,  Points  1-5 
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Figure  A-16:  Profile  V-4L  -  August  7th,  Set  B,  Points  6-9 
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Fluorescence  (ppb)  q  |  Fluorescence  (ppb) 


i  A-17:  Profile  H-5L  -  August  7th,  Set  B,  Points  11-14 


Figure  A-18:  Profile  H-IM  -  August  2nd,  Set  A,  Points  5-10 
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Fluorescence  (ppb)  ^  |  Fluorescence  (ppb) 


i  A-19:  Profile  V-IM  -  August  2nd,  Set  A,  Points  10-14 


Figure  A-20:  Profile  H-2M  -  August  2nd,  Set  B,  Points  3-6 
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Fluo 


Figure  A-21:  Profile  V-2M  -  August  2nd,  Set  B,  Points  9-13 
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Fluorescence  (ppb) 


Figure  A-23:  ProfQe  V-3M  -  August  5th,  Set  A,  Points  5-8 
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Fluorescence  (ppb)  q  |  Fluorescen 


A-25:  Profile  V-4M  -  August  5th,  Set  B,  Points  5-9 
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Figure  A-27;  Profile  H-IH  -  July  31st,  Set  B,  Points  2-6 


Figure  A-28:  ProfUe  H-2H  -  July  31st,  Set  C,  Points  1-5 
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Figure  A-29:  Profile  V-IH  -  July  31,  Set  C,  Points  5-9 
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Figure  A-30:  Profile  H-IZ  -  July  24tli,  Set  C,  Points  2-5 
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Fluorescence  (ppb) 
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Figure  A-31:  Profile  H-2Z  -  July  24th,  Set  C,  Points  6-10 


Figure  A-32:  Profile  V-IZ  -  August  8th,  Set  A,  Points  10-12 
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Fluorescence  (ppb) 


Figure  A-33:  Profile  V-2Z  -  August  Sth,  Set  A,  Points  12-14 


Figure  A-34:  Profile  H-3Z  -  August  Sth,  Set  A,  Points  15-19 
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Figure  A-35:  Profile  V-3Z  -  August  8th,  Set  A,  Points  20-23 
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Appendix  B  -  Daily  Data  Summaries 


Data  is  presented  for  each  day  of  field  study.  Setups  on  each  day  are  designated  as  Run 
A,  Run  B,  and  so  on.  Data  is  also  included  on  current  and  water  depth  measurements. 


July  24*  Data 


Run  C 

Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

3:57:00 

100.00 

55.00 

20.60 

49.40 

661.00 

3.77 

2 

4:04:00 

80.00 

55.00 

38.40 

62.20 

118.00 

11.20 

3 

4:06:00 

90.00 

55.00 

51.70 

76.80 

91.00 

15.80 

4 

4:08:20 

100.00 

55.00 

16.40 

30.20 

91.00 

6.20 

5 

4:10:20 

110.00 

55.00 

0.16 

0.13 

86.00 

0.03 

6 

4:12:45 

110.00 

45.00 

2.49 

6.20 

76.00 

1.39 

7 

4:14:00 

100.00 

45.00 

20.10 

29.50 

91.00 

6.07 

8 

4:16:00 

90.00 

45.00 

134.00 

93.90 

91.00 

19.30 

9 

4:18:30 

80.00 

45.00 

142.00 

71.70 

91.00 

14.70 

10 

4:20:30 

70.00 

45.00 

49.50 

73.50 

91.00 

15.10 

11 

4:22:30 

80.00 

35.00 

18.60 

33.00 

91.00 

6.79 

12 

4:24:30 

90.00 

35.00 

10.10 

18.90 

91.00 

3.88 

13 

4:26:30 

100.00 

35.00 

16.70 

17.50 

91.00 

3.59 

14 

4:28:30 

110.00 

35.00 

19.50 

26.10 

241.00 

3.30 

Currents 

Point 

Time 

Current 

Sigma  Readings 

Error 

Depth 

1 

3:45 

7.01 

1.17  3 

1.97 

27" 

2 

•  4:00 

7.61 

0.96  5 

0.91 

3 

4:30 

11.53 

0.26  6 

0.21 

CM 

CM 

SaiDDle  Time 

=  3:00 

Sample  Time  =  4:00 

Sample  Time  =  4:30 

t  (seconds)  d(cm) 

u  (cm/sec) 

t  (seconds)  d  (cm) 

u  (cm/sec) 

t  (seconds)  d(cm) 

u  (cm/sec) 

16.71 

100 

5.98 

10.92 

100  9.16 

8.7 

100 

11.49 

14.80 

100 

6.76 

15.14 

100  6.61 

8.43 

100 

11.86 

12.06 

100 

8.29 

12.91 

100  7.75 

8.65 

100 

11.56 

13.57 

100  7.37 

8.8 

100 

11.36 

13.97 

100  7.16 

8.52 

100 

11.74 

8.97 

100 

11.15 

Average 

7.01 

Average 

7.61 

Average 

11.53 

Sigma 

1.17 

Sigma 

0.96 

Sigma 

0.26 

Points 

3 

Points 

5 

Points 

6 

97 


July  31®*  Data 


Run  B 

Point 

Time 

H.  Coord. 

V.  Coord. 

F  (ppb) 

sigma  F 

#  of  pts. 

Error 

1 

9:44:00 

110.00 

30.00 

164.98 

49.68 

30 

17.78 

2 

9:45:00 

100.00 

30.00 

133.53 

16.29 

30 

5.83 

3 

9:46:00 

90.00 

30.00 

219.86 

5.04 

30 

1.80 

4 

9:47:00 

70.00 

30.00 

186.29 

8.37 

30 

3.00 

5 

9:48:00 

50.00 

30.00 

40.70 

15.35 

30 

5.49 

6 

9:49:00 

40.00 

30.00 

0.65 

0.17 

30 

0.06 

7 

9:50:00 

50.00 

30.00 

0.35 

0.08 

30 

0.03 

8 

9:51:00 

80.00 

30.00 

0.57 

0.38 

30 

0.14 

Run  C 

Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

10:02:00 

80.00 

35.00 

0.26 

0.26 

90 

0.05 

2 

10:04:00 

100.00 

35.00 

56.63 

28.50 

90 

5.89 

3 

10:06:00 

110.00 

35.00 

59.78 

9.31 

90 

1.92 

4 

10:08:00 

120.00 

35.00 

2.05 

1.42 

90 

0.29 

5 

10:10:00 

105.00 

35.00 

66.74 

15.83 

30 

5.66 

6 

10:11:00 

105.00 

25.00 

260.41 

13.55 

90 

2.80 

7 

10:13:00 

105.00 

15.00 

92.24 

18.36 

90 

3.79 

8 

10:15:00 

105.00 

5.00 

1.89 

0.91 

90 

0.19 

9 

10:17:00 

105.00 

45.00 

1.37 

1.48 

260 

0.18 

10 

10:22:00 

105.00 

30.00 

57.56 

52.37 

270 

6.25 

11 

10:27:00 

125.00 

30.00 

45.92 

25.13 

270 

3.00 

12 

10:32:00 

145.00 

30.00 

3.12 

12.71 

315 

1.40 

13 

10:37:45 

110.00 

20.00 

74.45 

17.43 

360 

1.80 

14 

10:44:15 

110.00 

40.00 

7.07 

5.29 

467 

0.48 

Currents 

Point  Time 

Current 

Sigma 

Readings 

Error 

Depth 

1  9:25 

2.32 

0.16 

7 

0.12 

2  9:53 

3.18 

0.22 

6 

0.18 

33" 

3  10:20 

4.09 

0.34 

5 

0.32 

30" 

4  10:45 

6.37 

0.89 

7 

0.65 

27" 

98 


July  31®*  Data  (cont.) 


Sample  Time  -  9:25 


t  (sec) 

d  (inches) 

u  (cm/sec) 

19.00 

16.00 

2.14 

16.40 

16.00 

2.48 

6.02 

6.00 

2.53 

6.27 

6.00 

2.43 

7.08 

6.00 

2.15 

6.65 

6.00 

2.29 

6.85 

6.00 

2.22 

Average 

Sigma 

Points 

Sample  Time  =10:20 

2.32 

0.16 

7 

t  (sec) 

d  (inches) 

u  (cm/sec) 

6.51 

12.00 

4.68 

7.80 

12.00 

3.91 

7.81 

12.00 

3.90 

7.91 

12.00 

3.85 

7.40 

12.00 

4.12 

Average  4.09 

Sigma  0.34 

Points  5 


Sample  Time  =  9:53 


(sec) 

d  (inches) 

u  (cm/sec) 

4.56 

6.00 

3.34 

5.50 

6.00 

2.77 

4.58 

6.00 

3.33 

4.57 

6.00 

3.33 

4.91 

6.00 

3.10 

4.75 

6.00 

3.21 

Average 

Sigma 

Points 

Sample  Time  =  10:45 

3.18 

0.22 

6 

t  (sec) 

d  (inches) 

u  (cm/sec) 

5.47 

12.00 

5.57 

5.07 

12.00 

6.01 

6.03 

12.00 

5.05 

4.57 

12.00 

6.67 

4.58 

12.00 

6.66 

4.37 

12.00 

6.97 

3.97 

12.00 

7.68 

Average 

6.37 

Sigma 

0.89 

Points 

7 
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August  2“**  Data 


Run  A 


Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

11:35:00 

68.00 

40.00 

140.25 

32.89 

90.00 

6.80 

2 

11:37:00 

55.00 

40.00 

122.52 

46.33 

90.00 

9.57 

3 

11:39:00 

45.00 

40.00 

9.42 

9.70 

90.00 

2.00 

4 

11:41:00 

35.00 

40.00 

8.29 

4.98 

30.00 

1.78 

5 

11:42:00 

35.00 

40.00 

56.82 

6.92 

30.00 

2.47 

6 

11:43:00 

25.00 

40.00 

4.59 

0.90 

30.00 

0.32 

7 

11:44:00 

45.00 

40.00 

157.77 

25.72 

30.00 

9.20 

8 

11:45:00 

55.00 

40.00 

195.60 

15.58 

30.00 

5.57 

9 

11:46:00 

65.00 

40.00 

4.96 

1.50 

30.00 

0.54 

10 

11:47:00 

50.00 

40.00 

200.57 

3.20 

30.00 

1.15 

11 

11:48:00 

50.00 

30.00 

93.47 

5.68 

30.00 

2.03 

12 

11:49:00 

50.00 

20.00  ■ 

6.93 

1.55 

30.00 

0.55 

13 

11:50:00 

50.00 

50.00 

75.25 

1.88 

30.00 

0.67 

14 

11:51:00 

50.00 

60.00 

2.25 

0.23 

30.00 

0.08 

15 

11:52:00 

50.00 

40.00 

95.15 

35.27 

30.00 

12.62 

Run  B 

Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

12:01:00 

85.00 

40.00 

2.30 

8.32 

210.00 

1.12 

2 

12:05:00 

70.00 

40.00 

5.66 

6.02 

150.00 

0.96 

3 

12:08:00 

60.00 

40.00 

91.92 

23.80 

150.00 

3.81 

4 

12:11:00 

50.00 

40.00 

33.36 

12.87 

150.00 

2.06 

5 

12:14:00 

40.00 

40.00 

15.55 

21.09 

210.00 

2.85 

6 

12:18:00 

70.00 

40.00 

12.32 

14.91 

210.00 

2.02 

7 

12:22:00 

50.00 

40.00 

44.36 

19.57 

90.00 

4.04 

8 

12:24:00 

50.00 

30.00 

30.12 

35.47 

30.00 

12.69 

9 

12:25:00 

70.00 

30.00 

67.73 

7.98 

90.00 

1.65 

10 

12:27:00 

70.00 

20.00 

2.94 

1.38 

90.00 

0.28 

11 

12:29:00 

70.00 

40.00 

31.62 

18.49 

90.00 

3.82 

.  12 

12:31:00 

70.00 

50.00 

1.70 

3.45 

390.00 

0.34 

13 

12:38:00 

70.00 

35.00 

46.06 

18.84 

210.00 

2.55 

14 

12:42:30 

70.00 

30.00 

13.23 

8.86 

200.00 

1.23 

Currents 

Point  Time 

Current 

Sigma 

Readings 

Error 

Depth 

1  11:20 

1.34 

0.18 

8 

0.12 

27" 

2  11:55 

3.61 

0.41 

15 

0.19 

3  12:15 

5.04 

0.42 

8 

0.28 

25.5" 

4  12:35 

3.35 

0.21 

10 

0.12 

24" 

100 


August  2“**  Data  (cont.) 


Sample  Time  =  11:20 
t  (seconds)  d  (inches)  u  (cm/sec) 


13.7 

6 

1.11 

13.5 

6 

1.13 

12.5 

6 

1.22 

11.5 

6 

1.33 

10.5 

6 

1.45 

10.65 

6 

1.43 

10.45 

6 

1.46 

9.5 

6 

1.60 

Average  1 .34 

Sigma  0.18 

Points  8 


Sample  Time  =12:15 


t  (seconds) 

d  (inches) 

u  (cm/sec) 

4.7 

9 

4.86 

4.8 

9 

4.76 

4.9 

9 

4.67 

5.1 

9 

4.48 

4 

9 

5.72 

4.2 

9 

5.44 

4.5 

9 

5.08 

4.3 

9 

5.32 

Average  5.04 

Sigma  0.42 

Points  8 


Sampie  Time  =  1 1 :55 
t  (seconds)  d  (inches)  u  (cm/sec) 


4.27 

6 

3.57 

4.29 

6 

3.55 

5.25 

6 

2.90 

4.23 

6 

3.60 

4.28 

6 

3.56 

4.76 

6 

3.20 

3.67 

6 

4.15 

5 

9 

4.57 

6.64 

9 

3.44 

7.23 

9 

3.16 

4.26 

6 

3.58 

4.26  ' 

6 

3.58 

3.72 

6 

4.10 

4.09 

6 

3.73 

4.34 

6 

3.51 

Average 

3.61 

Sigma 

0.41 

Points 

15 

Sample  Time  =  12:35 


t  (seconds)  d  (inches) 

u  (cm/sec) 

11.3 

15 

3.37 

11.3 

15 

3.37 

12.1 

15 

3.15 

12.2 

15 

3.12 

11.8 

15 

3.23 

11.2 

15 

3.40 

11.4 

15 

3.34 

11.1 

15 

3.43 

9.9 

15 

3.85 

11.8 

15 

3.23 

Average 

3.35 

Sigma 

0.21 

Points 

10 
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August  5***  Data 


Run  A 


Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

2:25:00 

100.00 

40.00 

37.55 

13.70 

30.00 

4.90 

2 

2:26:00 

110.00 

40.00 

296.83 

7.15 

30.00 

2.56 

3 

2:27:00 

120.00 

40.00 

55.20 

6.81 

30.00 

2.44 

4 

2:28:00 

130.00 

40.00 

0.60 

0.52 

30.00 

0.19 

5 

2:29:00 

110.00 

40.00 

306.76 

12.08 

30.00 

4.32 

6 

2:30:30 

110.00 

30.00 

155.28 

33.47 

30.00 

11.98 

7 

2:31:30 

110.00 

20.00 

0.10 

0.03 

30.00 

0.01 

8 

2:32:30 

110.00 

50.00 

8.98 

4.85 

30.00 

1.74 

Run  B 

Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

2:40:00 

110.00 

50.00 

22.40 

22.60 

90.00 

4.68 

2 

2:42:00 

120.00 

50.00 

57.70 

14.60 

90.00 

3.01 

3 

2:44:00 

130.00 

50.00 

122.00 

24.50 

90.00 

5.07 

4 

2:46:00 

140.00 

50.00 

10.50 

3.91 

90.00 

0.81 

5 

2:48:00 

133.00 

50.00 

77.00 

31.70 

90.00 

6.56 

6 

2:50:00 

133.00 

40.00 

72.10 

9.26 

90.00 

1.91 

7 

2:52:00 

133.00 

30.00 

4.18 

1.74 

90.00 

0.36 

8 

2:54:00 

133.00 

45.00 

127.00 

27.60 

90.00 

5.71 

9 

2:56:00 

133.00 

60.00 

5.88 

2.51 

270.00 

0.30 

Run  C 

Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

3:02:00 

130.00 

45.00 

67.40 

28.00 

210.00 

3.79 

2 

3:06:00 

140.00 

45.00 

29.90 

14.80 

210.00 

2.00 

3 

3:10:00 

120.00 

45.00 

0.98 

1.20 

210.00 

0.16 

4 

3:14:00 

130.00 

45.00 

61.60 

30.10 

210.00 

4.07 

5 

3:19:00 

130.00 

55.00 

21.30 

3.32 

210.00 

0.45 

6 

3:23:00 

130.00 

65.00 

0.99 

0.18 

90.00 

0.04 

7 

3:25:00 

130.00 

60.00 

4.53 

1.20 

90.00 

0.25 

8 

3:27:00 

130.00 

50.00 

73.00 

8.17 

90.00 

1.69 

9 

3:29:00 

130.00 

40.00 

55.70 

32.70 

90.00 

6.77 

Currents 

Point 

Time 

Current 

Sigma 

Readings 

Error 

Depth 

1 

2:11 

1.12 

0.02 

2 

0.09 

20" 

2 

2:35 

2.55 

0.27 

6 

0.22 

19" 

3 

3:00 

4.23 

0.14 

7 

0.1 

18" 

4 

3:15 

4.36 

0.49 

5 

0.47 

5 

3:30 

4.34 

0.22 

7 

0.16 

14" 
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August  5*  Data  (cont.) 


Sample  Time  =  2:11 

Sample  Time  =  2:35 

t  (seconds)  d  (inches) 

u  (cm/sec) 

t  (seconds)  d  (inches) 

u  (cm/sec) 

26.8  12 

1.14 

12.9 

12 

2.36 

27.4  12 

1.11 

12.1 

12 

2.52 

12.5 

12 

2.44 

12.9 

12 

2.36 

11.9 

12 

2.56 

9.9 

12 

3.08 

Average 

1.12 

Average 

2.55 

Sigma 

0.02 

Sigma 

0.27 

Points 

2 

Points 

6 

Sample  Time  =  3:00  Sample  Time  =  3:15 


t  (seconds) 

d  (inches) 

u  (cm/sec) 

t  (seconds)  d  (inches) 

u  (cm/sec) 

7 

12 

4.35 

6.83  12 

4.46 

7.6 

12 

4.01 

7.3  12 

4.18 

7.14 

12 

4.27 

7.73  12 

3.94 

7.07 

12 

4.31 

7.55  12 

4.04 

7.48 

12 

4.07 

5.9  12 

5.17 

7.25 

12 

4.20 

7 

12 

4.35 

Average 

4.23 

Average 

4.36 

Sigma 

0.14 

Sigma 

0.49 

Points 

7 

Points 

5 

Sample  Time  =  3:30 


t  (seconds)  d  (inches) 

u  (cm/sec) 

7.38 

12 

4.13 

7.43 

12 

4.10 

7.24 

12 

4.21 

7.17 

12 

4.25 

6.88 

12 

4.43 

6.51 

12 

4.68 

6.68 

12 

4.56 

Average 

4.34 

Sigma 

0.22 

Points 

7 
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August  6***  Data 


Run  A 

Point  Time 

1  2:27:00 

2  2:35:30 

3  2:36:30 

4  2:37:30 

5  2:39:00 

Run  B 

Point  Time 

1  3:29:00 

2  3:30:00 

3  3:31:00 

4  3:32:00 

5  3:33:00 

6  3:34:00 

7  3:35:00 

8  3:36:00 

9  3:37:00 

Run  C 


H.  Coord.  V.  Coord. 
120.00  45.00 

170.00  45.00 

180.00  45.00 

190.00  45.00 

190.00  35.00 


H.  Coord.  V.  Coord. 
150.00  50.00 

140.00  50.00 

130.00  50.00 

120.00  50.00 

110.00  50.00 

125.00  50.00 

130.00  50.00 

130.00  40.00 

130.00  60.00 


F  (ppb)  sigma  F 
71.77  91.60 

92.43  26.69 

83.44  35.62 

80.14  42.57 

42.72  70.24 


F  (ppb)  sigma  F 
0.05  0.00 

0.05  0.01 

3.12  1.15 

91.35  8.48 

3.29  2.52 

25.00  18.06 

51.99  18.27 

19.05  11.00 

3.58  2.23 


#  of  pts.  Error 
480.00  8.20 

30.00  9.55 

30.00  12.74 

30.00  15.23 

390.00  6.97 


#  of  pts.  Error 
30.00  0.00 

30.00  0.00 

30.00  0.41 

30.00  3.03 

30.00  0.90 

30.00  6.46 

30.00  6.54 

30.00  3.93 

30.00  0.80 


Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

3:45:00 

135.00 

60.00 

48.61 

9.04 

30.00 

3.23 

2 

3:46:00 

135.00 

50.00 

125.10 

27.15 

30.00 

9.72 

3 

3:47:00 

135.00 

40.00 

54.83 

6.42 

30.00 

2.30 

4 

3:48:00 

135.00 

30.00 

0.86 

0.24 

30.00 

0.09 

5 

3:49:00 

135.00 

50.00 

201.57 

9.16 

30.00 

3.28 

6 

3:50:00 

130.00 

50.00 

77.94 

22.23 

30.00 

7.95 

7 

3:51:00 

120.00 

50.00 

3.59 

2.95 

30.00 

1.06 

8 

3:52:00 

140.00 

50.00 

85.23 

30.51 

30.00 

10.92 

9 

3:53:00 

145.00 

50.00 

32.64 

16.80 

30.00 

6.01 

10 

3:54:00 

150.00 

50.00 

3.69 

1.42 

30.00 

0.51 

11 

3:55:00 

135.00 

50.00 

116.92 

61.09 

30.00 

21.86 

Currents 

Point 

Time 

Current 

Sigma 

Readings 

Error 

Depth 

1 

2:30 

1.35 

0.13 

5 

0.12 

20.5" 

(@  2:40) 

2 

3:20 

2.02 

0.27 

5 

0.26 

3 

3:40 

6.20 

0.70 

7 

0.51 

17.5" 

4 

4:00 

8.04 

0.62 

10 

0.36 

16.25" 
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August  6***  Data  (cont.) 


Sample  Time  =  2:30  Sample  Time  =  3:20 


t  (seconds)  d  (inches) 

u  (cm/sec) 

t  (seconds) 

d  (inches) 

u  (cm/sec) 

23.3 

12 

1.31 

12.4 

12 

2.46 

25.8 

12 

1.18 

17.1 

12 

1.78 

21.8 

12 

1.40 

16.8 

12 

1.81 

19.8 

12 

1.54 

15 

12 

2.03 

22.8 

12 

1.34 

15.2 

12 

2.01 

Average  1.35 

Sigma  0.13 

Points  5 


Average  2.02 

Sigma  0.27 

Points  5 


Sample  Time  =  3:40 


Sample  Time  =  4:00 


t  (seconds) 

d  (inches) 

u  (cm/sec) 

t  (seconds) 

d  (inches) 

u  (cm/sec) 

5.6 

12 

5.44 

4.15 

12 

7.34 

5.8 

12 

5.26 

3.75 

12 

8.13 

4.9 

12 

6.22 

4.08 

12 

7.47 

5.1 

12 

5.98 

3.92 

12 

7.78 

4.7 

12 

6.49 

4.22 

12 

7.22 

4.4 

12 

6.93 

3.44 

12 

8.86 

4.3 

12 

7.09 

3.65 

12 

8.35 

3.85 

12 

7.92 

3.37 

12 

9.04 

3.66 

12 

8.33 

Average 

6.20 

Average 

8.04 

Sigma 

0.70 

Sigma 

0.62 

Points 

7 

Points 

10 
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August  7***  Data 


Run  A 

Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

4:25:00 

140.00 

30.00 

23.73 

13.67 

30.00 

4.89 

2 

4:26:00 

130.00 

30.00 

94.18 

24.20 

30.00 

8.66 

3 

4:27:00 

120.00 

30.00 

18.00 

10.76 

30.00 

3.85 

4 

4:28:00 

110.00 

30.00 

0.27 

0.20 

30.00 

0.07 

5 

4:29:00 

110.00 

50.00 

0.16 

0.05 

30.00 

0.02 

6 

4:30:00 

120.00 

50.00 

42.14 

7.60 

30.00 

2.72 

7 

4:31:00 

125.00 

50.00 

118.77 

23.14 

30.00 

8.28 

8 

4:32:00 

130.00 

50.00 

191.97 

37.97 

30.00 

13.59 

9 

4:33:00 

140.00 

50.00 

60.77 

13.32 

30.00 

4.77 

10 

4:34:00 

150.00 

50.00 

2.90 

0.68 

30.00 

0.24 

11 

4:35:00 

125.00 

50.00 

177.92 

20.05 

60.00 

5.07 

12 

4:36:00 

125.00 

60.00 

74.82 

9.86 

30.00 

3.53 

13 

4:37:00 

125.00 

40.00 

73.42 

26.06 

60.00 

6.59 

14 

4:38:00 

125.00 

30.00 

5.45 

1.26 

30.00 

0.45 

15 

4:39:00 

125.00 

50.00 

28.79 

39.67 

30.00 

14.20 

Run  B 

Point 

Time 

H.  Coord. 

V.  Coord. 

F(ppb) 

sigma  F 

#  of  pts. 

Error 

1 

4:44:00 

125.00 

50.00 

1.90 

6.28 

90.00 

1.30 

2 

4:46:00 

135.00 

50.00 

61.33 

38.36 

90.00 

7.93 

3 

4:48:00 

140.00 

50.00 

143.29 

35.79 

90.00 

7.40 

4 

4:50:00 

145.00 

50.00 

168.03 

44.06 

90.00 

9.10 

5 

4:52:00 

155.00 

50.00 

36.82 

44.76 

90.00 

9.25 

6 

4:54:00 

140.00 

50.00 

146.59 

40.56 

90.00 

8.38 

7 

4:56:00 

140.00 

40.00 

32.36 

8.79 

90.00 

1.82 

8 

4:58:00 

140.00 

60.00 

125.13 

39.96 

90.00 

8.26 

9 

5:00:00 

140.00 

65.00 

23.23 

11.74 

90.00 

2.42 

10 

5:02:00 

140.00 

55.00 

71.57 

62.15 

90.00 

12.84 

11 

5:12:00 

150.00 

55.00 

50.76 

22.68 

210.00 

3.07 

12 

5:16:00 

130.00 

55.00 

1.33 

2.57 

210.00 

0.35 

13 

5:20:00 

120.00 

55.00 

0.30 

0.88 

210.00 

0.12 

14  5:24:00 

Currents 

140.00 

55.00 

151.78 

28.88 

210.00 

3.91 

Point 

1 

Time 

4:04 

4:15 

Current 

2.55 

Sigma 

0.31 

Readings 

12 

Error 

0.16 

Depth 

18.5" 

2 

4:42 

7.21 

0.66 

7 

0.48 

17" 

3 

5:02 

7.66 

0.73 

10 

0.42 

14.5" 

4 

5:17 

7.06 

0.62 

8 

0.42 

13" 
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August  7***  Data  (cont.) 


Sample  Time  =  4:15  Sample  Time  =  4:42 


t  (seconds)  d  (inches) 

u  (cm/sec) 

t  (seconds)  d  (inches) 

u  (cm/sec) 

13.44 

12 

2.27 

4.8 

12 

6.35 

9.8 

12 

3.11 

4.6 

12 

6.63 

10.23 

12 

2.98 

4.37 

12 

6.97 

10.11 

12 

3.01 

3.8 

12 

8.02 

12.8 

12 

2.38 

3.94 

12 

7.74 

12.04 

12 

2.53 

4.43 

12 

6.88 

12.04 

12 

2.53 

3.86 

12 

7.90 

12.74 

12 

2.39 

13.38 

12 

2.28 

12.84 

12 

2.37 

12.89 

12 

2.36 

12.99 

12 

2.35 

Average 

2.55 

Average 

7.21 

Sigma 

0.31 

Sigma 

0.66 

Points 

12 

Points 

7 

Sample  Time  = 

5:02 

Sample  Time 

=  5:17 

t  (seconds)  d  (inches) 

u  (cm/sec) 

t  (seconds)  d  (inches) 

u  (cm/sec) 

3.69 

12 

8.26 

4.08 

12 

7.47 

3.62 

12 

8.42 

3.79 

12 

8.04 

3.45 

12 

8.83 

4.55 

12 

6.70 

4.12 

12 

7.40 

4.31 

12 

7.07 

4.44 

12 

6.86 

3.97 

12 

7.68 

4.34 

12 

7.02 

4.55 

12 

6.70 

3.85 

12 

7.92 

4.75 

12 

6.42 

4.37 

12 

6.97 

4.79 

12 

6.36 

4.47 

12 

6.82 

3.79 

12 

8.04 

Average 

7.66 

Average 

7.06 

Sigma 

0.73 

Sigma 

0.62 

Points 

10 

Points 

8 
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August  8***  Data 


Run  A 

Point 

Time 

H.  Coord. 

V.  Coord. 

F  (ppb) 

sigma  F 

#  of  pts. 

Error 

1 

5:06:00 

85.00 

45.00 

11.04 

15.66 

630.00 

1.22 

2 

5:17:00 

95.00 

45.00 

6.69 

11.04 

330.00 

1.19 

3 

5:23:00 

85.00 

45.00 

5.30 

6.20 

90.00 

1.28 

4 

5:25:00 

95.00 

45.00 

6.98 

9.45 

90.00 

1.95 

5 

5:27:00 

105.00 

45.00 

8.28 

7.87 

90.00 

1.63 

6 

5:29:00 

115.00 

45.00 

3.65 

4.73 

270.00 

0.56 

7 

5:34:00 

125.00 

45.00 

4.91 

7.64 

150.00 

1.22 

8 

5:37:00 

135.00 

45.00 

1.20 

1.75 

90.00 

0.36 

9 

5:39:00 

145.00 

45.00 

0.11 

0.23 

90.00 

0.05 

10 

5:41:00 

125.00 

45.00 

8.43 

7.51 

90.00 

1.55 

11 

5:43:00 

125.00 

40.00 

54.35 

24.42 

90.00 

5.05 

12 

5:45:00 

125.00 

36.00 

148.61 

47.27 

90.00 

9.77 

13 

5:47:00 

125.00 

45.00 

19.78 

15.15 

90.00 

3.13 

14 

5:49:00 

125.00 

40.00 

37.05 

34.57 

150.00 

5.53 

15 

5:52:00 

115.00 

40.00 

36.80 

26.58 

.  90.00 

5.49 

16 

5:54:00 

105.00 

40.00 

11.69 

28.73 

90.00 

5.94 

17 

5:56:00 

125.00 

40.00 

88.42 

40.13 

90.00 

8.29 

18 

5:58:00 

135.00 

40.00 

19.16 

22.61 

90.00 

4.67 

19 

6:00:00 

145.00 

40.00 

1.58 

1.69 

90.00 

0.35 

20 

6:02:00 

125.00 

40.00 

1j22.83 

^51.69 

38.45 

90.00 

7.94 

21 

6:04:00 

125.00 

45.00 

22.07 

90.00 

4.56 

22 

6:06:00 

125.00 

50.00 

16.67 

14.29 

90.00 

2.95 

23 

6:08:00 

125.00 

55.00 

4.09 

4.00 

90.00 

0.83 

Currents 

Point  Time 

Current 

Sigma 

Readings 

Error 

Depth 

1 

5:00 

3.16 

0.33 

9 

0.20 

16" 

2 

5:20 

6.14 

0.47 

11 

0.26 

15.5" 

3 

5:35 

8.06 

1.33 

9 

0.82 

14" 

4 

6:07 

13.26 

0.86 

7 

0.63 

11.25" 
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August  8***  Data  (cont.) 


Sample  Time  =  5:00 
t  (seconds)  d  (inches)  u  (cm/sec) 


11.7 

12 

2.61 

10.76 

12 

2.83 

10.11 

12 

3.01 

9.93 

12 

3.07 

9.83 

12 

3.10 

9.3 

12 

3.28 

8.65 

12 

3.52 

8.87 

12 

3.44 

8.42 

12 

3.62 

Average  3.16 

Sigma  0.33 

Points  9 


Sample  Time  =  5:35 
t  (seconds)  d  (inches)  u  (cm/sec) 


4.37 

12 

6.97 

4.83 

12 

6.31 

4.64 

12 

6.57 

3.09 

12 

9.86 

3.63 

12 

8.40 

3.81 

12 

8.00 

3.55 

12 

8.59 

3.92 

12 

7.78 

3.04 

12 

10.03 

Average  8.06 

Sigma  1.33 

Points  9 


Sample  Time  =  5:20 
t  (seconds)  d  (inches)  u  (cm/sec) 


5.55 

12 

5.49 

4.84 

12 

6.30 

5.39 

12 

5.65 

5.17 

12 

5.90 

5.4 

12 

5.64 

4.91 

12 

6.21 

4.99 

12 

6.11 

4.92 

12 

6.20 

4.88 

12 

6.25 

4.46 

12 

6.83 

4.36 

12 

6.99 

Average 

6.14 

Sigma 

0.47 

Points 

11 

Sample  Time  -  6:07 


t  (seconds)  d  (inches) 

u  (cm/sec) 

2.39 

12 

12.75 

2.11 

12 

14.45 

2.51 

12 

12.14 

2.35 

12 

12.97 

2.16 

12 

14.11 

2.42 

12 

12.60 

2.21 

12 

13.79 

Average  13.26 

Sigma  0.86 

Points  7 
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